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ABSTRACT 


The  average  radiation  flux  and  flux  divergence  in  the  atmosphere  during  March  is  computed  as  a function 
of  the  chief  absorbents  of  terrestrial  and  solar  radiation,  and  the  physical  parameters  of  the  atmosphere  such 
as  pressure,  temperature  and  cloudiness.  The  infrared  flux  for  water  vapor  is  computed  graphically  from 
Elsasser  radiation  charts.  The  absorption  of  insolation  by  water  vapor  is  computed  with  the  aid  of  an 
empirical  formula.  The  net  radiation  flux  divergence  is  then  computed  by  methods  of  finite  differences. 

It  is  shown  that  the  entire  troposphere  constitutes  an  energy  source  insofar  as  radiative  processes  are 
concerned,  and  that  at  all  latitudes,  the  largest  radiational  heat  loss  is  in  the  middle  troposphere.  It  is  also 
shown  that  the  effect  of  cloudiness  is  to  concentrate  this  maximum  radiational  cooling  to  layers  at  about  the 
height  of  middle  cloudiness  (3  to  4 km).  The  net  heat  loss  at  this  level  is  found  to  be  about  1.5  to  2.0°C 
day-1.  It  is  suggested  that  the  major  contribution  to  heating  the  troposphere,  to  offset  the  loss  due  to  radia- 
tion, is  the  release  of  latent  heat  by  condensation. 
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PREFACE 


The  major  portion  of  this  paper  was  originally  submitted  to  the  Graduate  School  of  Arts  and  Sciences, 
l\ew  York  University,  as  a doctoral  thesis.  The  results  reported  represent  the  first  completed  phase  of 
research  on  the  heat  balance  project  under  sponsorship  of  the  Geophysics  Research  Directorate  of  the  Air 
Force  Cambridge  Research  Center  (Contract  No.  AF  19(122)-165). 

It  was  the  original  purpose  of  the  research  to  investigate  the  relationship  between  the  latitudinal  dis- 
tribution of  radiative  excess  and  deficiency  and  various  proposed  models  of  the  mean  meridional  circulation 
of  the  atmosphere  during  a transitional  month  (March).  This  purpose  has  since  been  revised  to  afford  a 
more  complete  study  of  the  radiative  heat  loss  in  the  atmosphere.  At  the  present  time,  the  research  is 
being  continued  along  two  lines.  The  first  is  to  extend  the  present  study  to  all  seasons  of  the  year.  The 
second  is  to  include  other  factors  besides  the  radiative  components  which  contribute  to  the  atmospheric 
heat  budget.  It  is  hoped  that  the  completion  of  this  next  phase  of  research  will  in  part  supplement  our 
understanding  of  the  driving  mechanism  of  the  circulations  of  the  atmosphere. 

I wish  to  thank  Professor  II.  A.  Panofskv,  formerly  of  New  York  University  and  now  at  Pennsylvania 
State  College,  for  his  valuable  aid  and  constant  encouragement  during  all  phases  of  the  reported  research, 
and  Professor  J.  Spar  of  New  York  University  for  his  painstaking  review  of  the  original  manuscript. 
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THE  DISTRIBUTION  OF  RADIATIONAL  TEMPERATURE 
CHANGE  IN  THE  NORTHERN  HEMISPHERE  DURING  MARCH* 

1.  INTRODUCTION 

All  energy  of  atmospheric  motion  ig  ultimately  derived  from  incoming  solar  radiation.  Methods  for  the 
transformation  of  insolational  energy  to  that  of  mean  motion  of  the  atmosphere  are  discussed  by  Bjerknes 
et  al.  (1933)  and  Rosshy  (1941)  among  others. 

Since  the  mean  temperature  of  the  earth  and  its  atmosphere  changes  very  little  from  year  to  year, 
incoming  solar  radiation  must  be  balanced  by  radiation  going  out  to  space,  if  both  radiation  currents  are 
integrated  over  a period  of  a year  over  both  hemispheres.  At  any  particular  time  and  place,  however,  there 
need  not  be,  and  probably  is  not,  a radiation  balance.  The  transformation  of  energy  of  radiation  to  that  of 
atmospheric  motion  is  a consequence  of  the  inequality  of  incoming  and  outgoing  radiation  currents  as  func- 
tions of  time  and  space. 

The  radiation  unbalance  at  the  top  of  the  atmosphere  results  from  an  excess  of  incoming  over  outgoing 
radiation  in  the  equatorial  regions,  and  a deficiency  in  the  polar  regions.  Because  of  the  large  difference 
between  the  temperature  of  the  sun  and  that  of  the  earth  and  its  atmosphere,  the  spectral  distributions  of  the 
energy  of  solar  (short-wave)  and  terrestrial  (long-wave)  radiation  are  almost  completely  mutually  exclusive. 
As  a result,  the  atmosphere  is  semitransparent  to  solar  radiation  but  almost  opaque  to  terrestrial  radiation. 

Only  a small  portion  of  the  energy  of  the  incoming  solar  beam  is  used  directly  to  heat  the  atmosphere. 
The  rest  is  either  reflected  hack  to  space,  or  is  used  to  heat  the  earth’s  surface  and  evaporate  moisture  from 
the  large  bodies  of  water.  The  energy  which  reaches  the  earth’s  surface  is  then  transferred  to  the  atmosphere 
by  conduction  and  eddy  convection  in  the  form  of  sensible  heat  and  in  the  form  of  latent  heat.  These 
constitute  the  major  heat  sources  for  the  atmosphere. 

It  was  shown  by  Mtfller  and  Miigge  (1931)  that  the  atmosphere,  at  least  through  the  troposphere,  con- 
tinually loses  heat  through  long-wave  radiation.  Since  the  absorption  of  insolation  by  the  atmosphere  is 
comparatively  small,  the  atmosphere  provides  a heat  source  insofar  as  radiative  processes  are  concerned. 
This  was  first  clearly  indicated  by  Moller  (1935).  Subsequent  investigations  (Elsasser  (1940  a),  Tanck 
(1940)),  while  limited  in  scope,  have  substantiated  this  fact. 

The  purpose  of  this  study  is  to  present  a unified  picture  of  the  latitudinal  and  vertical  distribution  of 
average  radiative  heat  loss  in  the  northern -hemisphere  troposphere.  The  study  covers  both  clear  skies  and 
conditions  of  average  cloudiness.  This  is  to  provide  an  estimate  of  the  importance  of  cloudiness  to  the 
problem  of  the  atmospheric  radiation  budget.  The  month  of  March  was  chosen  for  the  research  because 
the  radiation  budget  of  the  atmosphere  during  March  in  many  ways  resembles  the  average  annual  budget. 
It  is  hoped  that,  in  the  future,  similar  studies  can  be  made  for  the  other  months  of  the  year,  thus  providing 
a more  complete  description  of  the  role  played  by  the  atmosphere  as  a radiative  heat  source. 


Manuscript  received  for  publication  30  July  1952. 
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A few  words  might  be  included  at  this  point  concerning  the  methods  used  in  calculating  the  short-wave 
and  long-wave  radiation  fluxes  and  the  net  flux  divergence. 

The  depletion  of  the  solar  beam  has  been  studied  experimentally  by  Fowle  (1915),  Kimball  (1927)  and 
others,  and  theoretically  by  Yamamoto  and  Onishi  (1949  a,  1951).  Their  results  are  in  good  agreement 
with  those  of  Fowle  and  Kimball.  The  calculations  of  absorption  of  insolation  by  water  vapor  in  the  atmos- 
phere are  straightforward,  albeit  tedious. 

The  computations  of  long-wave  radiation  flux,  however,  are  complicated  by  the  fact  that  the  infrared 
absorption  spectrum  of  water  vapor  is  extremely  variable.  For  this  reason  most  previous  investigators  such 
as  Simpson  (1928),  and  others  have  had  to  make  simplifying  assumptions  concerning  the  absorption  coeffi- 
cients of  water  vapor  which,  although  certainly  valid  for  establishing  the  correct  order  of  magnitude  of  the 
net  outward  flux  at  the  top  of  the  atmosphere,  could  not  show  up  the  correct  detailed  distribution  of  the 
radiation  flux. 

A graphical  method  of  integration  of  the  basic  radiative  transfer  equations  was  first  developed  by  Mtlgge 
and  Midler  (1932  b),  thus  providing  for  a more  exact  computation  of  the  net  radiation  flux  in  the  atmosphere. 
With  increased  knowledge  of  the  water  vapor  absorption  spectrum,  came  more  sophisticated  methods  of 
computing  radiation  flux.  Elsasser  (1940  b)  developed  a new  radiation  chart  by  devising  a method  of 
expressing  the  extremely  variable  water  vapor  absorption  coefficients  by  an  analytic  function  which  was 
both  simple  and  accurate  enough  to  be  of  practical  value. 

The  difficulties  encountered  in  interpreting  and  applying  spectroscopic  theory  and  data  to  water  vapor 
absorption  in  the  atmosphere  have  led  to  a more  direct  experimental  approach  to  the  problem  of  infrared 
radiation  in  the  atmosphere.  Emissivity  measurements  in  the  laboratory  and  the  atmosphere  made  by 
1 . A.  Brooks  (1941),  Elsasser  (1941),  Robinson  (1947,  1950)  and  others,  indicate  certain  empirical  relation- 
-hip«  between  emissivity  and  total  path  length  of  the  eminent.  These  relationships  were  used  by  D.  L. 
Brooks  (1Q18)  and  Robinson  (1950)  to  construct  radiation  charts  which  are  quite  similar  to  those  developed 
by  Mugce  and  M oiler,  and  by  Elsasser.  As  a matter  of  fact,  Elsasser  (1942)  and  D.  L.  Brooks  have  both 
shown  the  equivalence  of  the  radiation  charts  constructed  on  the  basis  of  spectroscopic  theory  and  observa- 
tion- to  those  constructed  by  using  emissivity  measurements  primarily.  The  present  computations  of 
infrared  radiation  flux  make  use  of  the  Elsasser  radiation  chart,  because  it  is  felt  that,  despite  some  of  the 
incomplete  spectroscopic  resolution  of  the  infrared  spectrum,  it  provides  the  more  general  basis  for  the 
theoretical  and  empirical  solution  to  the  transfer  problem  in  the  atmosphere. 

There  are  indications,  based  on  a study  by  Yamamoto  and  Onishi  (1949b),  that  minor  revisions  to  the 
K l»a«ser  Radiation  Chart  are  needed.  The  effect  of  these  revisions  are  estimated  by  Yamamoto  (1950). 
Other  recent  comments  and  criticisms  of  F.lsasser’s  solutions  to  the  transfer  problem  are  discussed  in  Sec- 
tion 2 of  this  report. 

The  net  radiative  flux  divergence  computed  in  this  study  is  derived  from  the  smoothed  vertical  dis- 
tribution of  the  radiation  flux.  The  values  of  the  temperature  changes  were  then  calculated  from  the 
flux  divergence  between  discrete  levels. 

Methods  have  been  suggested  by  Bruinenberg  (1946),  D.  L.  Brooks  (1950),  Hales  (1950)  and  Cowling 
(l'*50  b)  for  computing  point  radiative  temperature  changes  in  the  atmosphere  from  the  distribution  of  the 
flux  divergence.  These  methods,  however,  are  of  the  nature  of  mathematical  refinements,  consisting  essen- 
tially of  devices  for  smoothing  discrete  sets  of  data.  In  the  case  of  average  flux  distribution,  they  are  prob- 
ably no  more  accurate  than  the  more  primitive  methods  used  in  the  present  calculations. 


2.  THE  THEORY  OF  INFRARED  RADIATIVE  FLUX  IN  THE  ATMOSPHERE 


S.L  THE  RADIATIVE  TRANSFER  EQUATION 

The  basic  equation  for  downward  directed  monochromatic  radiative  transfer  through  an  absorbing  alab 
in  the  atmosphere  is  given  by  the  Schwaraschild  equation  (e.g.,  Brunt  (1939)) 

dl,  - — k,  sec  [(I.  - Ib,)  du,  (1) 

where  J,  is  the  intensity  of  the  incoming  radiation  beam  at  the  frequency  r;  k,  is  the  absorption  coefficient 
per  unit  density  per  unit  length  of  the  absorbent;  f is  the  zenith  angle  of  the  incoming  beam;  Ib,  is  the 
monochromatic  black -body  radiation  emerging  from  the  slab  and  du,  the  optical  depth  of  the  slab,  is  defined 
by  du  » pdz,  where  p is  the  density  of  the  absorbent  and  z is  height  designated  as  positive  upward. 

The  solution  to  the  transfer  equation  for  the  total  infrared  flux  was  given  by  Elaasaer  (1940  b)  as 

F - /*£%'***  (2) 


where  T is  the  temperature, /»,,  the  specific  black-body  flux,  r,  is  a function  representing  the  transmission  of 
radiative  flux  and  /„  a generalized  absorption  coefficient,  is  the  value  of  k,  averaged  over  a small  interval  of 
the  infrared  spectrum. 


One  of  the  chief  difficulties  in  evaluating  Eq.  (2)  is  that,  while  in  theory  it  should  be  possible  to  determine 
the  form  of  the  transmission  function  77,  in  practice  this  is  no  simple  task. 


The  absorption  spectra  of  atmospheric  gases  such  as  water  vapor,  carbon  dioxide  and  ozone  are,  in 
general,  band  spectra  where  each  of  the  absorption  bands  is  composed  of  a large  number  of  narrow  absorption 
lines.  The  relationship  between  the  absorption  coefficient  and  the  frequency  for  these  spectral  lines  is  of  the 
form  shown  schematically  in  Fig.  1 (e.g.,  Dennison  (1928)).  In  Fig.  1,  r0  represents  the  frequency  at  the 
center  of  the  line,  and  the  half  width  a is  defined  as  the  frequency  interval  for  which  the  absorption  coefficient 
at  the  center  of  the  line  has  been  reduced  to  one  half  of  its  value.  The  total  intensity  of  the  spectral  line 
is  given  by 


S 


k(p)  dp. 


(3) 


Ladenburg  and  Reiche  (1913)  have  shown  that  in  the  case  of  a single  absorption  line  (or  widely  spaced 
lines),  the  absorption  takes  place  in  the  central  portion  of  the  line  for  small  optical  depths  and  is  proportional 
to  the  product  of  the  optical  depth  and  the  total  line  intensity.  As  the  optical  depth  becomes  larger,  how- 
ever, tb*  absorption  becomes  complete  at  the  center  of  the  line  and  further  absorption  takes  place  in  the  line 
wings,  where  (p  — »»)>■  0 snd  is  proportional  to  the  square  root  of  the  optical  depth,  the  total  line  intensity 
and  the  half  width  of  the  absorption  line.  The  fractional  transmission  is  related  to  the  absorption  by 
r m (1  — A)  where  A is  the  fractional  absorption  in  the  frequency  interval  Ar. 

The  above  analysis  has  been  extended  by  Elsasser  (1938,  1942)  to  the  case  of  absorption  in  a band 
spectrum  where  the  band  is  composed  of  an  infinite  train  of  equally  spaced  lines,  each  of  the  same  half  width 
and  total  intensity.  The  transmission  function  for  a beam  of  radiation  is  then  given  in  terms  of  the  error 
function  of  argument  y/lu/2,  where  the  generalized  absorption  coefficient  l is  defined  as  l «■  (2 waS)/tP. 
Here  a and  S are  the  half  width  and  total  intensity  of  each  absorption  line,  respectively,  and  d represents  the 
distance  between  lines.  In  the  case  of  absorption  through  a horizontal  layer  in  the  atmosphere,  it  is  shown 
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Fig.  1.  Shape  of  a spectral 
line  (schematic). 


that  the  transmission  function  for  the  radiation  flux  becomes 

“ [1  — erf  (y/lu/2)g  (hi/2)],  (4) 

where  g(lu/ 2)  is  an  algebraic  function  which  can  be  evaluated  numerically.  If  Elsasser’s  model  is  reduced 
to  a single  absorption  line,  Eq..(4)  can  be  reduced  to  the  solution  of  Ladenburg  and  Reiche. 

More  recently,  Matossi,  Mayer  and  Rauscher  (1949)  have  discussed  the  total  absorption  for  overlapping 
lines  of  equal  half  widths  but  different  total  intensities  and  distances  between  the  line  centers.  This  problem 
is  solved  by  an  approximate  method,  and  here  again  it  is  shown  that,  for  a single  absorption  line,  the  approxi- 
mate solution  for  overlapping  lines  reduces  to  the  exact  solution  given  by  Ladenburg  and  Reiche.  It  is  also 
shown  that  in  the  case  of  large  values  of  the  path  length,  the  square-root  law  of  Ladenburg  and  Reiche  is 
certainly  valid  for  overlapping  lines  as  far  as  the  dependence  on  the  path  length  is  concerned. 

These  results  have  been  confirmed  experimentally  by  Chapman  (1950)  using  a low  dispersion  prism 
spectrometer  on  the  water  vapor  bands  in  the  near  infrared.  The  observations  indicated  that  the  relation- 
ship between  the  fractional  absorption  in  these  bands  and  the  total  water  vapor  at  constant  pressure  is  of  the 
form  given  by  Elsasser.  For  a large  range  of  values  of  the  water  vapor  length,  the  relationship  takes  the  form 
of  the  square-root  absorption  law. 

Cowling  (1950  a)  has  recently  derived  a semiempirical  expression  for  the  absorption  of  infrared  radiation 
by  water  vapor.  The  expression  has  the  form  of  a simple  logarithmic  law  and  it  has  been  pointed  out  by 
Cowling  (1950  b)  that  for  calculations  within  the  troposphere,  the  results  are  almost  identical  to  those 
of  Elsasser. 

Finally,  mention  should  be  made  of  the  recent  extension  of  Elsasser’s  work  by  Kaplan  (1952  a,  1952  b). 
Here  again,  the  derived  transmission  functions  give  almost  the  same  results  for  radiative  flux  computations 
in  the  lower  atmosphere  as  those  discussed  above. 
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Since  water  vapor  is  the  primary  absorbent  of  infrared  radiation  in  the  atmosphere,  the  solution  to  the 
radiative  transfer  equation  depends  mainly  on  the  determination  of  a general  absorption  law  for  the  water 
vapor  spectrum.  The  total  flux  can  then  be  obtained  by  graphical  integration  as  was  first  suggested  by 
Miigge  and  Moller  (1932  a),  or  by  numerical  integration  as  was  suggested  by  Lonnquist  (1950).  A radia- 
tion chart  for  the  calculation  of  the  total  radiation  flux  has  also  been  developed  by  Elsasser  (1940  b). 

It  should  be  noted  that  the  development  for  the  total  flux  is  independent  of  the  nature  of  the  absorbing 
medium  and  can  be  applied  to  any  absorbent  in  the  atmosphere.  Suitable  transmission  functions,  as  dis- 
cussed by  Craig  (1949)  and  Yamamoto  (1950),  can  be  used  in  Eq.  (2)  for  carbon  dioxide  and  ozone,  since 
their  infrared  absorption  spectra  are  somewhat  simpler  than  that  of  water  vapor.  At  atmospheric  pressures, 
these  transmission  functions  would  be  reduced  nearly  to  the  functions  for  exponential  absorption  given  by 
Elsasser  (1942). 

Although  the  problem  of  finding  the  exact  transmission  function  (or  functions)  for  the  complicated  infra- 
red spectrum  of  water  vapor  has  not  yet  been  solved  completely,  Elsasser’s  results  seem  to  be  the  most  prac- 
tical for  meteorological  studies  of  the  lower  atmosphere.  These  results,  in  the  form  of  the  radiation  chart, 
have  been  used  as  the  basis  for  most  of  the  infrared  flux  calculations  discussed  below.  The  details  of  radia- 
tion flux  computations  using  the  Elsasser  chart  are  discussed  by  Elsasser  (1942)  and  Charney  (1945). 


2.2.  PRESSURE  AND  TEMPERATURE  CORRECTIONS 

The  expression  for  the  total  flux  given  in  Eq.  (2)  was  derived  under  the  assumption  that  the  absorption 
and  emission  in  the  atmosphere  are  independent  of  pressure,  and  dependent  on  temperature  only  through 
the  black-body  laws.  Unfortunately,  this  is  not  the  case. 

As  was  shown  above,  for  most  values  of  optical  path  in  the  atmosphere,  the  total  absorption  in  a small 
wavelength  interval  is  dependent  on  the  square  root  of  both  the  optical  depth  and  the  half  width  of  the 
absorption  line.  Although  the  total  intensity  of  the  absorption  line  is  independent  of  pressure  and  tempera- 
ture, the  shape  of  the  line,  shown  in  Fig.  1,  is  affected  by  a number  of  factors.  For  our  analysis,  the  most 
important  of  these  is  the  broadening  of  the  line  as  a result  of  collisions  between  air  molecules  and  the  mole- 
cules of  the  radiating  gas  in  the  atmosphere.  This  effect  is  called  pressure,  or  collision,  broadening  of  the 
spectral  lines.  As  a result  of  pressure  broadening,  the  line  becomes  flatter,  the  half  width  is  increased  and 
the  absorption  at  the  center  or  core  of  the  line  becomes  diminished,  although  there  is  an  increase  of  absorp- 
tion in  the  line  wings.  For  a single  absorption  line,  it  can  be  shown  from  kinetic  theory  that  the  half  width 
is  directly  proportional  to  the  pressure  and  inversely  proportional  to  the  square  root  of  the  temperature 
(e.g..  Brunt  1939). 

The  pressure  effect  in  the  line  wings  is  decreased  for  overlapping  absorption  lines,  such  as  those  studied 
by  Elsasser.  Thus,  a smaller  pressure  effect  should  be  expected.  Elsasser  (1942)  points  out  that  numerous 
investigators  have  shown  that  for  a large  range  of  the  infrared  spectrum,  in  the  cases  where  there  is  a suffi- 
cient amount  of  the  absorbent  (H20,  C02,  or  Oj),  so  that  the  square-root  law  of  absorption  holds,  the  total 
absorption  varies  as  the  fourth  root  of  the  air  pressure.  This  would  indicate  that  under  these  conditions  the 
half  width  is  proportional  to  the  square  root  of  the  pressure. 

This  square-root  pressure  correction  has  been  verified  by  Callendar  (1941)  and  Elsasser  (1942)  for  the 
15-M  band  of  carbon  dioxide  based  on  the  measurement  of  Hertz.  On  the  other  hand,  Moller  (1944)  has 
concluded  from  Wimmer’s  observation  of  the  4.3-ju  band  of  carbon  dioxide  that  insofar  as  absorption  is  a 
function  of  optical  path,  the  optical  path  for  carbon  dioxide  should  be  corrected  by  a factor  ip /pa)'11  where 
p0  is  standard  sea-level  pressure. 
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Strong  and  Watanabe  (1940),  and  Sumrnerfield  and  Strong  (1941  a,  1941  b)  have  also  indicated  that 
as  a result  of  overlapping  of  line  wings,  the  pressure  correction  for  infrared  absorption  of  ozone  and  water 
vapor  in  the  troposphere  is  of  the  form  .4  = bp'1  where  A represents  the  fractional  absorption  and  b is  a 
proportionality  constant.  This  result  is  implied  for  carbon  dioxide  as  well,  Chapman  (1950)  in  measuring 
the  fractional  absorption  in  the  1.35-p  and  1.85-ji  water  vaj>or  hands,  finds  a similar  pressure  correction. 
The  effect  of  pressure  on  the  6-m  band  of  water  vapor  has  been  investigated  by  \1atossi  and  Hauscher  (1949), 
who  similarly  find  that  the  total  absorption  is  very  nearly  proportional  to  the  fourth  root  of  pressure.  Their 
results,  however,  are  dependent  on  the  values  of  the  half  width  in  the  6-/j  hand  determined  by  Nielsen  (1941) 
of  0.52  cm'1.  This  half  width  is  somewhat  larger  than  that  given  by  Adel  (1947)  w ho  found  that  for  atmos- 
pheric water  vapor  lines  at  18.6  n and  15.99  n,  the  half  width  is  0.12  cm'1  and  0.1 1 cm  ',  respectively.  ’ Since 
the  smaller  the  half  width,  the  less  is  the  wing  effect  for  constant  optical  depth,  it  would  seem  that  there  is 
still  some  doubt  as  to  the  exact  correction  due  to  pressure  broadening. 

In  the  stratosphere  it  is  most  likely  that,  except  for  the  rotational  band  of  water  vapor  in  the  far  infra- 
red, the  individual  spectral  lines  become  well  separated  from  each  other  and,  as  suggested  In  Strong  and  Plass 
(1950),  a linear  pressure  correction  would  prevail.  The  amount  of  the  absorbents  in  the  stratosphere,  how- 
ever, is  quite  small,  and  the  difference  between  a linear  and  square-root  pressure  correction  would  he  impor- 
tant only  if  one  wanted  to  investigate  the  detailed  flux  and  flux  divergence  distribution  in  the  stratosphere. 
Since  this  study  is  confined  to  a model  principally  within  the  troposphere,  the  square-root  correction  is 
used  throughout/ 

The  radiation  flux  as  given  in  Eq.  (2)  is  a function  of  the  product  l u.  The  square-root  pressure  correc- 
tion can  therefore  be  expressed  as  a correction  to  the  optical  depth.  Thus  for  the  sake  of  convenience, 
we  write  

u*  = V p/p0  (5) 

where  u*  is  the  effective  optical  depth,  p is  the  average  pressure  in  the  layer  considered  and  p0  is  standard 
sea-level  pressure  (1013  mb). 

Very  little  is  known  about  the  relationship  between  line  shape  and  tcnqierature  in  the  case  of  hand 
spectra.  The  above  discussion,  however,  leads  to  the  conclusion  that  the  half  width  of  spectral  lines  in  the 
infrared  should  be  inversely  proportional  to  the  fourth  root  of  the  absolute  temperature  and  should  be  in- 
cluded in  the  correction  considered  in  Eq.  (5).  In  the  troposphere,  the  temperature  correction  would 
always  be  close  to  unity  and  it  has  been  omitted  in  the  calculations  below. 


3.  THE  PHYSICAL  MODEL  OF  THE  ATMOSPHERE  IN  THE  NORTHERN  HEMISPHERE 

DURING  MARCH 


It  has  been  shown  in  Section  2 that  the  distribution  of  radiation  flux  and  radiation  cooling  is  primarily 
a function  of  the  distribution  of  parameters  characterizing  the  physical  state  of  the  atmosphere.  The  model 
of  the  atmosphere  described  in  this  section  is,  for  the  most  part,  based  on  data  which  represent  conditions 
prevailing  during  average  March  in  the  northern  hemisphere.  Because  of  the  lack  of  data  in  the  Arctic, 


t More  recent  evidence  indicates  a swing  of  the  pendulum  back  to  a linear  pressure  correction.  It  can  be  shown,  however,  that 
this  would  not  seriously  affect  the  calculations  for  the  lower  and  middle  troposphere. 


17 


the  present  model  is  limited  to  the  region  0°  to  70s  N.  The  following  discussion  covers  the  distribution  with 
height  snd  latitude  of  atmospheric  pressure  and  temperature,  cloudiness  and  the  chief  absorbents  of  terres- 
trial  and  solar  radiation  such  as  water  vapor,  ozone  and  carbon  dioxide. 


3.1.  TEMPERATURE  AND  PRESSURE 
3.11.  Troposphere 

The  temperature  and  pressure  distributions  used  in  this  study  were  obtained  from  the  Northern  Hemi- 
sphere Upper  Level  Charts  for  March  published  by  the  U.  S.  Weather  Bureau  (1944).  The  data  on  these 
charts  sre  in  the  form  of  isopleths  of  temperature  and  pressure  for  the  entire  northern  hemisphere.  Tempera- 
ture and  pressure  values  were  averaged  from  these  charts  for  each  10  degrees  of  latitude  (0°  to  70°  N)  around 
complete  latitude  circles,  for  each  level  of  published  data.  From  these,  the  vertical  distributions  of  pres- 
sure and  temperature  were  calculated.  It  was  necessary  to  use  extrapolated  values  for  the  equatorial 
regions  since  the  isopleths  were  not  generally  drawn  south  of  about  10°  N.  For  this  reason,  the  tempera- 
tures for  the  0°  to  10°  N latitude  belt,  particularly  in  the  upper  levels,  were  adjusted  to  agree  with  the  data 
published  by  Wagner  (1931),  Ratner  (1946)  and  Hess  (1948). 


Fig.  2.  Temperature  distribu- 
tion in  the  troposphere  and 
overage  height  of  the  tropo- 
pause  during  March  (0°  to 
70°  N). 
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3.12.  The  Heigfu  of  the  Tropopauae 

The  tropopauae  was  considered  to  be  the  upper  bounding  "surface”  for  most  of  the  radiation  dux  com- 
putations. By  use  of  the  vertical  distribution  of  temperature  discussed  above,  it  is  possible  to  depict  an 
approximate  variation  with  latitude  of  the  height  of  the  tropopauae.  In  the  middle  latitudes,  the  tropo- 
pause  becomes  indistinct  and  in  fact  is  a relatively  deep  zone  of  transition  rather  than  a discontinuity  of  lapse 
rate.  The  average  heights  used  here  for  the  middle  latitudes  follow  from  the  definition  of  the  height  of 
the  tropopause  suggested  by  Flohn  (1947).  Flohn  defines  the  "height  of  the  tropopauae”  as  the  level  at 
which  the  average  lapse  rate  curves  above  and  below  the  "tropopause  layer”  intersect.  \ vertical  cross 
section  showing  the  distribution  of  temperature  and  the  variation  of  the  height  of  the  tropopause  during 
March  is  shown  in  Fig.  2. 

3.13.  Stratosphere 

The  radiation  flux  downward  from  the  stratosphere  comes  from  a layer  extending  from  the  tropopause 
to  a height  of  about  25  km  (e.g.,  Karandikar  (1946)).  If  we  assume  that  the  pressure  distribution  in  the 
stratosphere  is  such  that  the  pressure  is  reduced  by  ahout  one  half  for  each  5-km  elevation  (Penndorf  (1946  b), 
Crimminger  (1947)),  that  is, 

pk  = 

where  h is  in  km,  we  can  easily  calculate  the  approximate  average  pressure  for  this  radiating  layer.  Since 
the  tropopause  slopes  downward  toward  the  poles,  the  average  pressure  of  the  lower  stratosphere  (to  about 
25  km)  will  increase  with  increasing  latitude. 

The  number  of  calculations  and  observations  of  temperature  in  the  lower  stratosphere  has  steadily- 
increased  within  the  last  five  years.  (For  reference,  see  Best  et  al.  (1947),  Cox  et  al.  (1949),  Murgatroyd 
and  Clews  (1949)  and  Brasefield  (1950),  each  of  which  contains  ample  references  to  this  point.)  Although 
these  observations  and  calculations  are  not  consistent  in  either  time  or  place,  they  help  to  form  a pattern 
of  the  thermal  structure  of  the  lower  stratosphere  which,  for  the  most  part,  agrees  well  with  the  data  pre- 
sented by  Penndorf  (1946  a,  1946  b)  and  Chiplonkar  (1940).  During  all  seasons  the  lower  stratosphere  is 
either  isothermal  or  has  a slight  inversion  up  to  at  least  25  km.  The  average  temperature  in  the  stratosphere 
between  60°  and  70°  N is  about  25°C  higher  than  in  the  equatorial  belt.  The  radiation  flux  from  the  strato- 
sphere is  very  insensitive  to  changes  of  average  temperature  of  5°C  and/or  average  pressure  of  10  mb.  This 
is  primarily  because  of  the  low  concentrations  of  most  of  the  absorbents  at  these  high  levels.  (See  the 
sections  on  the  distribution  of  water  vapor,  ozone,  and  carbon  dioxide.) 


3.2.  THE  DISTRIBUTION  OF  CLOUDS 

Clouds  act  as  black-body  radiators  for  infrared  radiation,  and  varying  reflectors  for  visible  radiation. 
It  is  therefore  important  to  know  the  distribution  of  cloudiness  by  cloud  type  as  well  as  the  average  height 
and  thickness  for  each  type. 

The  most  authoritative  study  of  the  distribution  of  total  cloud  cover  for  both  northern  and  southern 
hemispheres  is  that  conducted  by  Brooks  (1927).  This  study  gives  the  average  total  cloud  cover  over  the 
land  and  over  the  oceans  for  each  10°  latitude  belt  for  each  month  of  the  year.  Unfortunately,  the  total 
cloud  cover  has  not  been  broken  down  into  groups  of  cloud  types. 
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For  the  present  paper,  two  basic  sources  were  used  to  determine  the  distribution  of  total  cloud  cover 
with  latitude  and  the  percentage  of  cloud  cover  by  each  cloud  type.  The  results  were  then  compared  with 
thoee  of  Brooks  as  a qualitative  check  on  the  distribution  of  total  cloudiness.  The  basic  sources  are: 

1.  Preliminary  Climatic  Atlas  of  the  World  (1943),  which  gives  the  total  cloudiness  by  months  in 
percent  of  total  sky,  for  about  200  selected  land  stations  over  the  northern  hemisphere. 

2.  Atlas  of  Climatic  Charts  of  the  Oceans  (1938),  which  gives  the  seasonal  frequency  of  total  cloud  cover 
and  of  cover  by  various  cloud  types  over  the  oceans. 

The  data  from  the  Climatic  Atlas  of  the  Oceans  were  averaged  over  the  winter  and  spring  seasons  to 
give  approximate  values  for  March.  Means  of  the  values  of  total  cloudiness  were  taken  for  5°  squares  of 
longitude  and  latitude  from  both  atlases,  and  from  these,  the  average  total  cloud  cover  for  each  10°  latitude 
belt  between  0°  and  70°  N was  computed.  The  following  table  shows  the  comparison  between  Brooks’ 
values  and  those  computed  in  the  present  study,  of  the  distribution  of  total  cloudiness,  in  percent,  for  March. 

If  the  total  cloudiness  is  weighted  by  the  relative  area  in  each  10°  belt,  the  average  cloud  cover  for  the 
region  0s  to  70°  N is  48.2  percent  according  to  Brooks’  data  and  50.0  percent  according  to  the  present  study. 
The  agreement  for  both  the  distribution  of  total  cloudiness  with  latitude  and  the  average  cloudiness  for  the 
region  0°  to  70°  N is  good,  particularly  in  view  of  the  fact  that  Brooks  used  data  for  only  about  one  half  the 
ocean  area  for  the  northern  hemisphere.  The  agreement  justifies  the  use  of  the  data  from  the  atlases 
mentioned  above  for  the  breakdown  of  the  total  cloudiness  into  percentage  amounts  of  each  cloud  type. 

Table  1.  The  dittr.lration  of  total  cloudiness  in  percentage  of  sky  cover  daring  March. 


Latitude  °N 

Brooks 

Present 

Study 

0°-10° 

Land,  L 

48 

55 

Oceans,  O 

50 

47 

L + O 

49 

49 

10°-20° 

L 

>9 

32 

O 

47 

45 

L+O 

39 

41 

20°-S0° 

L 

33 

34 

O 

50 

50 

L+O 

41 

44 

30° -40° 

L 

44 

46 

O 

55 

59 

L + O 

49 

54 

40°-50° 

L 

52 

53 

O 

68 

69 

L + O 

57 

60 

50°-608 

L 

55 

54 

O 

62 

63 

L+O 

57 

56 

60°-70° 

L 

52 

53 

O 

67 

— 

L+O 

54 

S3 

20 


The  percentage  of  each  cloud  type  was  calculated  from  the  data  given  in  the  Atlas  of  Climatic  Chart* 
with  the  aid  of  the  same  method  used  in  determining  total  cloud  cover.  In  this  calculation,  however,  the 
following  assumptions  were  made: 

1.  The  percentage  of  sky  covered  by  each  cloud  type  is  given  by  the  percentage  frequency  reported  for 
each  type  multiplied  by  the  percentage  of  total  cloud  cover. 

2.  There  is  no  overlap  of  clouds.  (That  is,  the  total  cloudiness  for  any  latitude  belt  is  the  same  aB  the 
total  of  the  observed  amounts  by  cloud  types  for  that  belt.) 

The  first  assumption  introduces  possible  sources  of  errors  for  the  following  reasons: 

1.  The  data  used  for  the  distribution  by  cloud  type  are  for  the  oceanic  regions  only  and  these  have 
been  used  to  describe  the  distributions  over  the  entire  latitudinal  belts.  This  is  less  serious  in  lower  than  in 
higher  latitudes,  since  the  ratio  of  ocean-to-land  surface  decreases  with  increasing  latitude.  The  computed 
distribution  of  low  cloudiness  is  affected  more  than  that  for  middle  or  high  cloudiness,  since  low  clouds  are 
more  influenced  by  local  terrain. 

2.  The  cloud  data  used  are  based  on  observations  at  1200  GMT.  Diurnal  variations  of  cloudiness, 
then,  are  not  reflected  in  the  averages  for  local  observations,  but  rather  in  the  cloudiness  distributions  about 
the  latitude  belt.  Also,  the  difficulty  of  observing  high  clouds  during  the  night  would  indicate  that  obser- 
vations over  the  Pacific  Ocean  (where  local  time  of  the  observations  is  near  2400)  do  not  show  the  true  cloudi- 
ness amounts  in  that  region. 

The  second  assumption  (no  overlap  of  clouds)  required  Blight  adjustments  of  the  observed  data  of  cloud 
type  amounts  since  values  for  total  cloudiness  were  computed  from  two  sets  of  complementary  data.  This 
modification,  however,  did  not  change  any  of  the  amounts  by  more  than  five  percent  of  their  values  and  in 
most  cases  the  change  was  much  less. 

For  convenience,  the  different  cloud  types  were  put  into  the  following  six  quasi-homogeneous  groups: 
(a)  cirrus,  cirrocumulus,  cirrostratus;  (b)  allostratus,  altocumulus;  (c)  stratus,  stratocumulus,  fracto- 
cumulus;  (d)  nimbostratus;  (e)  cumulus;  (f)  cumulonimbus.  The  following  table  shows  the  average  cloud 
amount  by  cloud  groups  for  10°  latitude  belts  for  March,  with  values  given  in  percentage  of  cloud  cover. 

It  is  seen  from  the  table  that  the  Ns  and  St  generally  increase  whereas  the  Cu  and  Cb  generally 
decrease  with  increasing  latitude. 


Tabic  2.  The  distribution  of  average  cloud  amounts  by  groups,  during  March.  (Values  are  given  in  percent  of  sky  covered 

by  each  cloud  group.) 


Latitude  °N 

Ci 

As 

Ns 

St 

Cu 

Cb 

Total 

0°-l0° 

12  1 

5 0 

6.3 

7.6 

13.5 

4 5 

49 

10c-20° 

10  3 

4 9 

4 1 

6 6 

11.2 

4.9 

41 

20  - 30° 

10  6 

6.2 

4.8 

8.4 

10.7 

3.3 

44 

30°-40° 

12  0 

7.1 

7.7 

14.1 

10  4 

2.7 

54 

40°-50° 

10  9 

7.3 

10  9 

18  2 

9.7 

3.0 

60 

50° -60° 

9 6 

6 3 

11  3 

16.9 

9.1 

2.8 

56 

60°-70° 

7 9 

6 0 

11.2 

17.0 

8.4 

2 5 

53 

21 


1 able  3.  Hit  dintri  billion  of  avorape  cloud  height*  \ty  groups  (lioightH  arf  tabulated  in  km). 


Latitude  °N 

Ci 

An 

Tops 

Ns 

Bases 

Si 

Tops 

Cu 

Bane* 

Tope 

Cb 

Bases 

0-10° 

9.8 

4.4 

4.4 

IS 

15 

2.3 

15 

5 0 

18 

IO°-20° 

10  4 

4 8 

4 8 

1 8 

18 

2 8 

18 

5.7 

2 3 

20°-30° 

10  4 

SO 

5 0 

2 0 

2 0 

3.0 

2.0 

6.0 

2.5 

30° -40° 

to  0 

4.7 

4.7 

18 

1.8 

2.8 

1.8 

5.4 

2.3 

40°-50° 

8.7 

4.0 

4 0 

1 3 

13 

2 3 

1.3 

4 7 

2 0 

S0°-60° 

7 6 

3 6 

3 6 

1.2 

12 

2 2 

12 

4 3 

1.7 

60°-70° 

7 2 

3 6 

3 6 

1 0 

1 0 

2.0 

1.0 

4 0 

15 

Very  few  studies  have  been  made  of  the  variation  of  cloud  height  with  type  and  latitude.  Siiring  (1931) 
has  given  an  estimate  of  the  variation  for  middle  and  high  clouds  based  on  observations  made  prior  to  1930. 
Tabular  summaries  are  presented  in  Shaw’s  Manual  oj  Meteorology  (1936)  and  in  the  Hann-Siiring  Lehrbuch 
der  Meteorologie  (1938),  of  average  summer  and  winter  heights  for  11  cloud  types  for  some  dozen  stations 
located  at  latitudes  ranging  from  6°  S to  70°  N.  Monthly  averages  of  the  heights  of  cloud  bases  are  gi  ven  in 
the  Airways  Meteorological  Atlas  published  by  the  U.  S.  Weather  Bureau  (1941)  for  about  55  stations  in  the 
United  States. 

Means  of  the  summer  and  winter  values  tabulated  by  Shaw  and  Hann-Siiring,  and  the  March  values 
tabulated  in  the  Airways  Meteorological  Atlas  were  plotted  for  each  cloud  group  as  functions  of  latitude. 
Smooth  curves  were  then  drawn  to  approximate  these  values,  and  the  curves  were  used  to  determine  the 
average  height  of  the  base  and  top  of  each  cloud  group  for  each  10°  latitude  belt  (from  0°  to  70°  N). 

It  is  obvious  that  a number  of  approximations  must  be  made  in  order  to  keep  the  cloud  data  in  a reason- 
able working  form  and  free  from  unimportant  and  cumbersome  detail.  For  the  present  study  it  is  assumed 
that: 

1.  The  various  cloud  types  can  be  put  into  six  quasi-homogeneous  groups  (as  discussed  above). 

2.  The  tops  and  bases  of  the  altostratus-altocumuius  groups  can  be  approximated  by  one  mean  height 
for  each  10°  latitude  belt. 

3.  Nimbostratus  cloud  tops  extend  to  the  altostratus  level. 

4.  Cumulus  clouds  have  depths  ranging  from  3.5  km  at  20°  to  30°  to  2.5  km  at  60°  to  70°. 

5.  Nimbostratus,  stratocumulus  and  cumulus  clouds  have  the  same  base  heights  in  any  latitude  belt. 

The  average  value  of  the  height  of  each  cloud  group  (Table  3)  is  a maximum  for  all  groups  at  latitudes 
20°  to  30°  N and  decreases  to  a minimum  at  60°  to  70°  N. 

3.3  THE  DISTRIBUTION  OF  WATER  VAPOR 
3.31.  Troposphere 

Water  vapor  is  the  chief  absorbent  of  infrared  radiation  in  the  atmosphere.  For  the  radiative  flux  calcu- 
lations, therefore,  it  is  necessary  that  the  water  vapor  distribution  with  both  height  and  latitude  be  deter- 
mined with  a fair  degree  of  accuracy.  Unfortunately,  most  observations  of  moisture  content  of  the  atmos- 
phere have  been  for  heights  below  5 or  6 km.  These  observations  would  suffice  for  a general  estimate  of  the 
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radiation  flux  because  of  the  controlling  nature  of  the  physical  parameters  of  the  lower  troposphere  on  the 
large-scale  radiation  budget.  The  detailed  water  vapor  distribution  throughout  the  troposphere  is  of  major 
importance,  however,  in  a study  of  radiative  temperature  changes  because  we  then  must  evaluate  the  diver- 
gence of  the  radiation  flux. 

Empirical  formulas  giving  the  vertical  distribution  of  vapor  pressure  with  height  have  been  discussed  by 
Shaw  (1936),  Hann-Siiring  (1938),  Brunt  (1940),  Yamamoto  (1949)  and  others.  These  formulas,  based  on 
observations  of  water  vapor  content  at  low  levels,  are  chiefly  of  the  form 

e*/e0  = f(h),  (6) 

where  e*  is  the  vapor  pressure  at  height  h,  e0  is  the  surface  vapor  pressure,  and  f{h)  is  an  exponential  function. 

During  the  course  of  the  present  study,  4-year  averaged  soundings  for  March  for  10  different  stations  in 
the  United  States  were  analyzed  in  an  attempt  to  fit  the  moisture  data  to  the  suggested  empirical  relation- 
ships discussed  above.  The  results  were  not  satisfactory.  Further  attempts  were  made  to  fit  the  data  by 
the  method  of  least  squares,  using  the  relationship  suggested  by  Hann-Siiring, 

e*a  = e0 10'<«*+“‘>,  (7) 

where  the  coefficients  a and  b were  determined  from  the  data.  It  was  found  that  the  coefficients  did  not  vary 
systematically  with  latitude,  and  in  particular,  the  coefficient  of  the  quadratic  term  varied  by  as  much  as 
2 orders  of  magnitude.  Thus,  even  by  varying  the  coefficients  for  each  latitude,  the  results,  particularly  in 
the  higher  levels,  were  not  sufficiently  accurate  for  use  in  the  radiation  flux  and  flux  divergence  computations. 
At  the  tropopause,  for  instance,  where  the  temperature  is  assumed  to  he  known,  the  relative  humidity 
derived  from  the  least-square  formula  varied  from  2 to  3 percent  to  over  300  percent  for  various  latitudes. 
An  analytic  expression  apparently  cannot  be  used  to  represent  accurately  the  vertical  distribution  of  water 
vapor,  particularly  in  the  upper  levels  of  the  troposphere. 

The  distribution  of  relative  humidity  in  the  atmosphere  shows  a much  smaller  vertical  variation  than 
does  that  of  vapor  pressure  or  specific  humidity.  A more  fruitful  approach  to  the  moisture  distribution 
problem  would  be  to  assume  a reasonable  vertical  distribution  of  relative  humidity  and  from  this,  since  the 
temperature  distribution  is  known,  to  derive  the  vertical  variation  of  water  vapor  directly. 

The  surface  distribution  of  relative  humidity,  as  a function  of  latitude,  is  given  by  Szava-Kovats  (1938) 
for  average  January  and  July.  The  data  represent  10-year  averages  of  observations  taken  at  about  6000  sta- 
tions spread  over  the  earth.  Since  Shands  (1949)  has  shown  that  the  distribution  of  water  vapor  in  the 
lower  atmosphere  in  March  is  better  approximated  by  winter  than  by  summer  conditions,  the  average  values 
of  surface  humidity  for  March  were  computed  by  weighting  the  January  values  3 times  as  much  as  those  in 
July.  Any  bias  in  the  data  because  of  the  diurnal  variation  of  relative  humidity  at  the  surface  was  elimi- 
nated by  the  process  of  averaging  the  data  over  complete  latitude  circles. 

Average  values  of  relative  humidity  for  March  for  levels  from  1 to  8 km  were  tabulated  from  the  data 
given  by  Ratner  (1946).  At  the  level  of  the  tropopause,  the  relative  humidity  was  assumed  to  be  near 
saturation.  This  assumption  is  supported  by  observational  evidence  published  by  Dobson  et  al.  (1946), 
Barrett  et  al.  (1949  a)  and  others.  The  values  tabulated  from  Ratner’s  data,  along  with  those  derived  for  the 
surface  and  assumed  for  the  tropopause,  were  plotted  as  functions  of  height  for  each  10°  latitude  belt.  A 
2-dimensional  field  of  relative  humidity  was  constructed  from  which  it  was  possible  to  deduce  the  relative 
humidity  distribution  in  the  upper  levels  of  the  troposphere.  The  vertical  and  latitudinal  distribution  of 
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Fig.  3.  Assumed  distribution 
of  relstive  humidity  (X)  in  the 
troposphere  daring  March. 


relative  humidity  (Fig.  3)  shows  an  axis  of  dry  air  at  about  9 to  10  lun  over  the  equator  to  about  6 km  at 
60°  to  70°  N.  Along  this  axis,  the  relative  humidity  changes  from  a minimum  of  about  25  percent  at  20°  N 
latitude  to  about  60  percent  at  60°  to  70°  N. 

The  distribution  of  specific  humidity  (and  total  water  vapor  path  length)  with  height  and  latitude 
through  the  troposphere  can  now  be  computed  since  the  temperature,  pressure  and  relative  humidity  are 
given.  As  was  indicated  in  the  preceding  paragraph,  since  the  variation  of  relative  humidity  in  the  atmos- 
phere is  much  less  than  the  variations  of  vapor  pressure  or  specific  humidity,  this  numerical  method  will 
give  more  accurate  results  for  the  moisture  content  distribution  than  would  the  approximation  formulas 
discussed  previously. 


The  total  water  vapor  content  between  any  two  levels  xt  and  xt  in  a vertical  column  of  unit  cross  section 
is  given  by 


u 


p.  ds  - (<j/g)(pt  ~ pi). 


(8) 


where  p»  is  the  density  of  the  water  vapor  in  the  layer,  p,  and  p\  are  the  pressures  at  the  bottom  and  top  of 
the  layer,  q is  the  mean  specific  humidity  in  the  layer,  and  g is  the  acceleration  of  gravity.  From  Eqs.  (5) 
and  (8),  the  effective  optical  depth  for  water  vapor  is 


«*  “ (q/g)^p/po  A/>, 


(9) 


where  p is  the  average  total  pressure  in  the  layer  A p,  and  p»  is  standard  sea -level  pressure. 


3.32.  Stratotphere 

As  a result  of  improved  instrumentation  in  recent  years,  there  have  been  a number  of  direct  observations 
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Table  4.  Summary  of  the  distribution  of  average  pressure,  temjHjrature  and  total  effective  water  vapor  content 

(to  the  trojHjpausc)  during  March. 
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of  the  humidity  distribution  in  the  lower  stratosphere  (e.g.,  Gliickauf  (1944),  (1945);  Brewer  et  al.  (1948); 
Barrett  et  al.  (1950)).  The  data  are  still  too  few,  however,  to  give  a reliable  picture  of  the  distribution  of 
the  amount  of  water  vapor  in  the  lower  stratosphere  with  latitude  or  season,  and  these  observations  can  at 
best  be  used  only  as  a guide  to  the  order  of  magnitude  of  the  stratospheric  moisture  content. 

In  his  studies  of  terrestrial  radiation,  Simpson  (1928)  assumed  that  the  stratosphere  contained  0.03  gm 
cm-*  of  precipitable  water.  On  the  other  hand.  Brunt  (1940)  and  Sheppard  (1947)  have  suggested  that  the 
total  water  vapor  above  the  tropopause  is  much  less,  that  is,  something  of  the  order  of  10_J  to  10“4  gm  cm"1. 
In  the  absence  of  more  complete  observational  evidence,  it  is  necessary  to  approximate  the  total  moisture 
content  using  reasonable  assumptions  based  on  our  knowledge  of  the  structure  of  the  lower  stratosphere. 

The  available  observations  show  that  the  relative  humidity  generally  increases  just  below  the  tropopause 
to  a maximum  at  the  tropopause  and  then  decreases  sharply,  remaining  less  than  10  percent  in  the  lower 
stratosphere.  Brewer  (1949)  points  out  that  this  sharp  decrease  in  relative  humidity  above  the  tropopause 
is  probably  due  to  the  circulation  in  the  lower  stratosphere,  lie  suggests  that  air  ascending  into  the  strato- 
sphere in  the  equatorial  regions  where  it  is  cold  and  dry  is  transported  to  temperate  and  polar  latitudes  where 
it  descends,  thus  accounting  for  a dry  stratosphere  at  all  latitudes. 

To  obtain  an  upper  limit  for  the  water  vapor  content  of  the  stratosphere,  it  is  reasonable  to  assume 
that  the  air  is  nearly  saturated  at  the  level  of  the  tropopause,  and  that  the  specific  humidity  in  the  lower 
stratosphere  is  constant  with  height.  Then,  as  suggested  by  Brunt  and  Kapur  (1938),  the  maximum  pre- 
cipitable water  in  the  stratosphere  would  depend  only  on  the  temperature  of  the  tropopause. 

The  total  mass  of  air  above  the  tropopause  per  unit  cross  section  is  given  by  M = />,/g,  where  p,  is  the 
pressure  at  the  tropopause.  The  specific  humidity  at  the  tropopause  is  given  approximately  by 

q,  » 0.62 (e./p,), 

where  e,  is  the  saturation  vapor  pressure  at  the  temperature  at  the  tropopause  level.  Then,  if  we  assume 
that  the  specific  humidity  is  constant  with  height  above  the  tropopause,  a maximum  value  for  the  total 
mass  of  water  vapor  in  a column  of  unit  cross  section  above  the  tropopause  (the  optical  depth)  is 

u = qt'M  S 0.62(e,/g).  (10) 

The  effective  optical  depth  for  the  stratosphere  is 

u*  = Vp/po  u, 

where  p is  the  average  total  pressure  in  the  lower  stratosphere.  In  the  computations  for  u *,  the  relative 
humidity  at  the  tropopause  level  has  been  assumed  as  90  percent  for  all  latitudes.  This  is  probably  an  over- 
estimate for  the  polar  regions,  and  it  would  have  the  effect  of  giving  slightly  too  large  downward  infrared 
radiation  in  latitudes  50°  to  60°  and  60°  to  70°  IN.  This  effect  has  been  investigated  by  Iloffer  (1951)  and 
found  to  be  quite  negligible. 

Table  4 given  below  shows,  in  summary  form,  the  distribution  of  pressure,  temperature  and  total  effec- 
tive optical  depth  (assuming  u*  = 0 at  the  height  of  the  tropopause)  used  in  the  present  study.  The  data, 
representing  average  conditions  for  March  for  each  10°  latitude  belt  (0°  to  70°  IN),  are  tabulated  for  each 
kilometer  up  to  the  tropopause,  as  well  as  for  the  average  downward  radiating  layer  in  the  lower  stratosphere. 
The  heigh  ts  of  the  bases  and  tops  of  the  various  cloud  groups  are  included  for  convenience  in  calculating  the 
radiation  flux  and  flux  divergence  in  the  troposphere. 
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Fig.  4.  Distribution  of  ozone 
density  (10~*  cm  Oj  km-1 
NTP)  in  tbe  stmosphere 
during  Msrcb. 
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3.4.  THE  DISTRIBUTION  OF  ATMOSPHERIC  OZONE 

The  ozone  concentration  in  the  atmosphere  has  rather  large  temporal  and  spatial  variations.  These 
variations  have  been  discussed  thoroughly  by  Gotz,  Meetham  and  Dobson  (1934),  Gotz  (1949)  and  Craig 
(1950).  It  is  well  established  that  there  is  a maximum  amount  of  ozone  in  the  northern  hemisphere  during 
March  and  that  the  total  ozone  content  increases  from  about  0.25  cm  (at  NTP)  at  the  equator  to  more  than 
0.30  cm  in  the  60°  to  70°  latitude  (at  NTP)  belt,  and  then  decreases  slightly.  The  average  values  for  the 
distribution  of  total  ozone  with  latitude  used  in  this  study  are  taken  from  the  study  of  Craig,  who  had  the 
benefit  of  a large  number  of  recent  observations  not  available  to  the  earlier  investigators. 
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The  vertical  distribution  of  ozone  has  been  measured  and  computed  by  means  of  Umkehr  effect,  balloon 
spectrograph  measurements,  photochemical  equilibrium  theory  and,  more  recently,  by  rocket  measurements. 
The  data  and  references  are  all  summarized  by  Craig.  These  measurements  and  computations  show  marked 
similarity  with  each  other  in  their  indicated  vertical  distributions,  and  it  has  been  indicated  by  Gotz  (1949) 
that  the  distribution  of  ozone  density  with  height  is  probably  a function  of  total  ozone  content.  The  value 
of  ozone  density  is  uniform,  and  quite  small  from  the  surface  layers  through  the  troposphere.  Above  the 
tropopause,  however,  the  ozone  density  increases  steadily  to  a maximum  at  about  25  km,  and  then  decreases 
above.  Measurements  made  from  studies  in  rocket  spectroscopy  show  no  significant  amount  of  ozone  above 
50  km  (Baum  et  al.  (1946)).  The  variation  shown  below  in  Fig.  4 is  based  primarily  on  the  distributions 
given  by  Gotz,  Meetham  and  Dobson  (1934)  and  Tonsberg  and  Olsen  (1944)  with  qualitative  changes  based 
on  recent  reports  of  Durant  and  Tousey  (1947)  and  Regener  (1950).  Included  in  Fig.  4 are  the  values  of  the 
distribution  of  total  ozone  content  with  latitude  used  in  the  present  computations. 


3.5.  THE  DISTRIBUTION  OF  CARBON  DIOXIDE 


Observations  of  carbon  dioxide  content  in  the  atmosphere,  which  have  been  made  since  at  least  1867 
(see  Callendar  (1940)),  do  not  indicate  any  systematic  variation  in  the  amount  of  carbon  dioxide  either  with 
season  or  latitude.  The  values  most  often  quoted  for  the  percentage  volume  of  carbon  dioxide  in  air  varies 
from  0.025  (Gliickauf  (1944))  to  0.030  (Callendar  (1940)).  On  the  basis  of  observations  made  in  England, 
Gliickauf  also  has  found  that  the  carbon  dioxide  content  of  the  free  air  is  constant,  at  least  from  4 to  10  km. 
Observations  of  the  composition  of  the  lower  stratosphere  (Paneth  (1939))  also  indicate  that  the  percentage 
volume  of  carbon  dioxide  in  air  is  constant  at  about  0.03  percent  at  least  up  to  30  km. 

If  we  assume  that  the  percentage  volume  of  carbon  dioxide  to  air  is  constant  with  height  and  latitude  at 
0.03  percent,  we  can  calculate  the  total  amount  of  carbon  dioxide  uc  in  a unit  vertical  column  from  the  surface 
to  the  top  of  the  atmosphere  from 

/QO 

pedz  = (mc/ m„)  X 3 X 10” 


where  mc/ma  is  the  ratio  of  the  molecular  weights  of  carbon  dioxide  to  air,  pe  and  p.  are  the  densities  of  carbon 
dioxide  and  air  in  a unit  volume  of  the  mixed  gas.  Equation  (11)  can  now  be  rewritten  as 


u,  = (me/ma)  X 


3 X 10- 

8 


4 rv c 

J. dp- 


(12) 


~ 0.465  gm  cm-*. 

If  the  total  carbon  dioxide  content  is  reduced  to  standard  pressure  and  temperature  this  becomes 


ue 


Pc 


I 


where  the  primes  indicate  values  at  standard  temperature  and  pressure  and 


Thus, 


Pc  — {pic/ m a) Pd  . 

0.465 

1.94  X 10”* 


(13) 

(14) 


240  cm. 


(15) 
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Table  S. 


The  vertical  distribution  of  total  effective  path  length  for  carbon  dioxide  in  the  atmosphere. 


p(mb) 

u,(gm  cm"*) 

u,'(cm  NTP) 

uc'*(cm  NTP) 

1000 

0.465 

240 

240 

900 

0.418 

215 

204 

800 

0.372 

192 

172 

700 

0.326 

168 

141 

600 

0.279 

144 

112 

500 

0.232 

120 

85 

400 

0.186 

96 

61 

300 

0.140 

72 

39 

200 

0.093 

48 

21 

100 

0.046 

24 

8 

75 

0.035 

18 

5 

50 

0.023 

12 
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The  absorption  by  carbon  dioxide,  as  in  the  case  of  water  vapor,  is  a function  of  the  total  pressure  and 
here  again  this  correction  can  be  applied  to  the  optical  path  for  carbon  dioxide.  Although  the  exact  nature 
of  the  pressure  dependency  is  still  somewhat  in  doubt  (see  Section  2),  a square-root  correction  has  been 
assumed  in  Table  5.  The  vertical  distribution  of  carbon  dioxide  and  effective  path  length  of  carbon  dioxide 
is  given  in  the  table. 


4.  CALCULATIONS  OF  INFRARED  FLUX  AND  INFRARED  COOLING 

IN  THE  TROPOSPHERE 

4.1.  COMPUTATION  OF  INFRARED  FLUX 

The  infrared  radiation  flux  computations  were  accomplished  by  means  of  Elsasser  radiation  charts 
(2nd  ed.).  For  the  sake  of  convenience,  these  computations  were  made  at  levels  (referred  to  as  reference 
levels)  determined  by  the  heights  chosen  as  representing  the  tops  and  bases  of  the  various  cloud  groups  in 
each  10°  latitude  belt.  This  was  done  to  avoid  unnecessary  duplication  in  the  computations  of  the  net  flux 
in  the  cases  of  clear  and  average  cloudy  skieB.  The  designation  of  these  reference  levels  is  given  in  Table  6. 

Radiation  charts  were  plotted  for  values  of  the  total  effective  optical  depth  against  the  vertical  distribu- 
tion of  temperature  for  each  of  the  seven  reference  levels  listed  above  for  each  10°  latitude  belt  from  0°  to 
70°  N.  The  net  radiation  flux  through  each  level  was  then  determined  by  planimetric  measurement  of  the 
areas  on  the  Elsasser  chart  representing  the  total  upward  and  the  total  downward  flux  through  the  reference 
level. 


Table  6.  Reference  levels  chosen  for  infrared  flux  computations. 


Level 

Remarks 

A 

ground 

B 

just  below  Stcu,  base  of  Ns  and  Cu 

C 

tops  of  Cu 

D 

just  above  As,  tops  of  Ns 

E 

tope  of  Cb 

F 

just  above  Ci 

G 

at  the  tropopause 

30 


35 


Fig.  6.  Distribution  of  net 
infrared  flux  (cal  cm'2  min*1) 
in  the  troposphere,  clear 
skies. 


Table  7.  The  distribution  of  net  infrared  flux  for  March,  clear  skies. 


0C 

’ to  J0°  N 

10" 

to  20°  N 

20° 

to  30°  N 

30" 

to  10"  N 

10 

to  50"  N 

50 

to  60  .\ 

60 

to 

70  N 

Level 

h 

F 

h 

F 

h 

F 

h 

F 

k 

F 

h 

F 

h 

h 

A 
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0 
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0 

0 

110 

0 
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0 

0 
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0 

0 
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0 

0 

l(S> 

0 

0 

no 

B 

l 

5 

0 
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1 8 

0 

159 

2 0 

0 172 

1.8 

0 

155 

1 3 

0 

133 

1 

2 

0 

130 

1 

0 

0 

122 

C 

2. 

3 

0 

160 

2 8 

0 

180 

3.0 

0 198 

2 8 

0 

176 

2 3 

0 

154 

2 

2 

0 

115 

2 

0 

0 

135 

D 

t 

4 

() 

201 

4 8 

0 

226 

5 3) 

0 239 

4 7 

0 

213 

4 0 

0 

190 

3 

6 

0 

173 

3. 

6 

0 

161 

E 

5 

0 

0 

213 

5 7 

0 

215 

6 0 

0 253 

5 4 

0 

227 

4 7 

0 

203 

4 

3 

0 

184 

4 

0 

0 

164 

F 

9 

8 

0 

297 

10  4 

0 

325 

10  4 

0 319 

10  0 

0. 

282 

8 7 

0 

247 

7 

6 

0 

215 

7 

2 

0 

186 

G(trop.) 

16 

0 

0 

354 

16  0 

0 

355 

16  0 

0 338 

13  0 

0 

2895 

no 

0 

255 

10 

0 

0 

226 

10 

0 

0 

196 

h is  the  height  of  the  reference  level  in  kilometers, 
f is  the  net  infrared  flux  in  cal  cm"*  min-1. 

4.11.  Clear  Skies 

The  net  flux  for  clear  skies,  calculated  for  seven  different  levels  for  each  latitude  belt,  is  given  in  Table  7. 
A curve  of  the  net  flux  as  a function  of  height  was  constructed  for  each  latitude  belt  (Figs.  5 (a)-(g))  and 
from  these  curves,  an  atmospheric  cross  section  (0°  to  70°  IS ) was  drawn  (Fig.  6)  showing  the  distribution 
with  height  and  latitude  of  the  net  infrared  radiation  flux  for  clear  skies. 

It  is  noted  from  Fig.  6 that  the  net  flux  is  everywhere  positive  — that  is,  directed  upward  — and  in- 
creases with  height.  The  net  flux  at  all  levels  is  largely  a function  of  the  temperature  at  the  earth’s  surface, 
since  some  25  to  30  percent  of  the  radiation  flux  through  the  tropopause  comes  directly  from  the  earth’s 
surface  and  this  percentage  is  even  larger  for  the  flux  through  the  lower  levels.  Other  conditions  being 
equal,  the  larger  the  water  vapor  content  of  the  atmosphere,  the  smaller  the  net  flux  at  all  levels,  because 
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there  is  then  a larger  downward  flux  to  balance  radiation  from  the  earth’s  surface.  Aloft,  an  increased  per- 
centage of  the  flux  comes  from  the  moist  layer  which  is  higher  and,  therefore,  colder  than  the  surface  of  the 
earth.  Thus,  the  combination  of  high  surface  temperature  and  low  atmospheric  humidity  results  in  a 
maximum  of  radiation  flux  at  about  15°  to  25°  N.  The  low  water  vapor  content  is  overbalanced  by  the  very 
low  surface  temperatures  to  give  minimum  flux  values  in  the  more  northerly  regions  (50°  to  70°  IV ). 

At  the  earth’s  surface,  with  clear  skies,  the  average  net  flux  in  units  of  cal  cm-2  min-1  is  about  0.110  with 
a slight  maximum  of  0.1 16  at  30°  to  40°  N.  At  a height  of  10  km  this  net  flux  increases  from  about  0.300  at 
the  equator  to  about  0.320,  then  decreases  to  about  0.195  at  higher  latitudes.  Since  the  tropopause  becomes 
lower  with  increasing  latitude,  the  net  upward  flux  through  the  tropopause  is  highest  in  equatorial  regions 
(0.350  at  16  km)  and  lowest  in  polar  regions  (about  0.200  at  10  km). 

4.12.  Modifications  Due  to  Cloudiness 

In  the  computations  described  above,  we  have  considered  the  net  infrared  flux  in  the  atmosphere  for 
conditions  of  clear  skies.  The  presence  of  clouds  in  varying  amounts  and  at  varied  heights  will,  of  course, 
alter  this  distribution  of  net  flux,  since  cloud  decks  radiate  and  absorb  as  black  bodies  in  the  infrared  region. 
Computations  for  clear  skies  must  then  be  modified  by  the  inclusion  of  the  black -body  radiation  upward  at 
the  temperature  at  the  tops  of  the  cloud  groups,  and  black -body  radiation  downward  at  the  temperature  at 
the  bases  of  these  cloud  groups.  If  one  is  now  to  compute  the  net  flux  through  any  level  aloft,  the  total 
upward  directed  flux  from  the  surface  must  be  modified  by  absorption  and  reradiation  (at  temperatures 
lower  than  at  the  earth’s  surface)  by  the  cloud  groups  below  the  level  considered.  Also,  the  downward  flux 
i6  increased  due  to  the  presence  of  black-bodv  surfaces  (cloud  decks)  above  the  reference  level  which  will 
absorb  infrared  sky  radiation  from  above  and  reradiate  this  energy  downward  through  the  reference  level. 

The  total  sky  is  only  partially  covered  by  each  cloud  group.  The  black-body  radiation  from  a cloud 
surface  through  any  level  in  the  atmosphere  is  therefore  decreased  by  an  amount  indicated  by  the  percentage 
of  skv  covered  by  that  cloud  surface.  Hence,  a cloud  group  covering  only  20  percent  of  the  sky  would  radiate 
one  fifth  of  the  black-bodv  radiation  of  a sky  completely  overcast  by  that  cloud  group.  The  effect  of  cloudi- 
ness is  to  reduce  the  net  flux  at  the  earth’s  surface  and  at  any  level  in  the  atmosphere. 

Calculations  of  net  flux  were  made  using  the  same  temperature-total  moisture  distributions  as  above, 
but  now  taking  into  account  the  average  cloudiness  for  each  latitude  belt.  The  results  of  these  calculations 


Table  8.  The  distribution  of  net  infrared  flux  for  March,  averaftp  cloudiness. 


0°  to  10°  N 

10°  to  20°  N 

20°  to  30°  N 

30°  to  40°  N 

40°  to  50 3 N 
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F 

h 

F 

h 

F 

h 

F 

h 

F 

h 

F 

\ 

0 

0 075 

0 

0 084 

0 
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0 

0 076 

0 

0 060 

0 

0 063 

0 

0 062 

B 

1 5 

0 096 

1 8 

0 117 

2.0 

0 125 

18 

0 098 

13 

0 073 

1.2 

0 071 

10 

0 066 

C 

2 3 

0 127 

2 8 

0 146 

3 0 

0 163 

2.8 

0 139 

2.3 

0 116 

2 2 

0 no 

2 0 

0 103 

D 

4 4 

0 172 

4 8 

0 199 

5 0 

0 210 

4 7 

0 185 

4 0 

0 168 

3 6 

0 151 

3 6 

0 143 
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5 0 

0 195 

5 7 

0 224 

6 0 

0 230 

5 4 

0 203 

4 7 

0 185 

4 3 

0 165 

4 0 

0 150 

F 

9 8 

0 271 

10  4 

0 296 

10  4 

0 282 

10  0 

0.249 

8 7 

0 222 

7 6 

0 197 

7 2 

0 173 

G(trop.) 

16  0 

0 325 

16  0 

0 324 

16  0 

0 304 

13  0 

0 256 

11  0 

0 226 

10  0 

0 207 

10  0 

0 184 

h is  the  height  of  the  reference  level  in  kilometers. 
F is  the  net  infrared  flux  in  cal  cm'*  min"1. 
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Fig.  8.  Distribution  of  net 
infrared  flux  (cal  cm*2  min"*) 
in  the  troposphere,  average 
cloudiness. 


are  given  in  Table  8.  A curve  of  the  net  flux  against  height  was  drawn  for  each  latitude  belt  (Figs.  7 (a)-(g)), 
and  from  these  curves,  values  were  obtained  showing  the  distribution  of  net  infrared  flux  (with  average  cloudi- 
ness) as  a function  of  height  and  latitude.  This  distribution  is  shown  in  Fig.  8. 

On  the  whole,  the  flux  distributions  for  clear  and  cloudy  skies  are  quite  similar.  The  net  flux  for  average 
cloudy  skies  is  everywhere  positive,  increases  with  height,  and  the  maximum  flux  at  each  level  is  again  at 
about  20°  to  25°  N . Although  the  net  flux  is  reduced  at  all  levels,  the  difference  in  the  net  flux  between  the 
ground  and  the  tropopause  is  increased  for  all  latitude  belts  (in  proportion  to  the  cloudiness  for  that  latitude 
belt).  For  instance,  it  is  found  that  the  net  flux  at  the  height  of  the  tropopause  is  reduced  by  about  0.020  to 
0.030  cal  cm'2  min'1  or  about  5 to  10  percent  of  its  value.  At  the  ground,  the  net  flux  is  reduced  by  about 
0.030  to  0.050  cal  cm'2  min'1  or  about  25  to  45  percent  of  its  value.  Thus,  the  presence  of  clouds  inhibits 
the  upward  flux  of  radiation  more  from  the  ground  than  from  levels  aloft. 

The  effectiveness  of  cloudiness  in  reducing  the  net  upward  flux  can  also  be  noted  when  comparing  Figs.  6 
and  8 according  to  latitude.  It  is  seen  that  at  the  ground,  the  smallest  reduction  is  found  at  latitudes 
10°  to  30°  IN  where,  not  only  is  the  total  average  cloudiness  a minimum  (41  to  44  percent  compared  with 
55  to  60  percent  at  latitudes  40°  to  60°  N)  but,  more  important,  the  amount  of  low  clouds  (stratus,  strato- 
cumulus  and  nimbostratus)  is  much  less  than  at  the  higher  latitudes  (22  to  24  percent  compared  with  35  to 
38  percent  in  the  more  northerly  regions).  The  effect  of  the  heights  of  the  various  cloud  groups  is  also  clearly 
seen  by  comparing  the  difference  of  the  net  flux  through  the  tropopause  with  clear  skies  against  that  where 
average  cloudiness  is  considered.  Here  the  reduction  due  to  cloudiness  is  at  a minimum  in  the  more  northerly 
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latitudes  where,  despite  the  larger  amount  of  total  cloudiness,  the  average  clouds  are  at  a lower  level  and 
therefore  radiate  upward  at  a higher  temperature  relative  to  the  surface  than  in  the  case  of  the  cloud  groups 
at  higher  levels  in  the  region  10°  to  30°  N. 

Another  very  significant  difference  between  the  two  distributions  is  the  fact  that,  for  clear  skies,  each 
graph  of  infrared  flux  against  height  (Figs.  5 (a)-(g))  has  a positive  and  almost  constant  slope  up  to  about 
5 to  8 km,  then  an  increasing  slope  to  the  height  of  the  tropopause.  In  the  case  of  average  cloudiness,  how- 
ever, the  analogous  graphs  (Figs.  7 (a )-(§))  have  positive  but  decreasing  slopes,  becoming  almost  horizontal 
for  heights  just  above  the  level  of  the  middle  clouds,  with  increasing  slopes  to  the  height  of  the  tropopause. 
This  difference  in  slopes  is  directly  related  to  the  difference  in  the  distribution  of  infrared  cooling  for  clear 
skies  and  for  average  cloudy  skies. 


4.2.  INFRARED  FLUX  DIVERGENCE 

It  has  long  been  recognized  that  for  both  clear  and  cloudy  skies,  the  earth’s  surface  and  its  troposphere 
continually  loses  heat  by  means  of  infrared  radiation  (Miigge  and  M oiler  1932  a,  Elsasser  1940  a).  The  heat 
loss  results  from  the  fact  that  most  of  the  radiation  leaving  the  earth’s  surface  is  absorbed  by  the  water  vapor 
and  clouds  above  and  then  re-emitted  both  upward  and  downward  at  temperatures  lower  than  at  the  ground. 
Since  the  amounts  of  both  water  vapor  and  cloudiness  decrease  with  height,  the  net  flux  for  a normal  tropo- 
sphere (temperature  decreasing  with  height)  is  directed  upward  and  increases  with  height  for  all  latitudes. 
Strong  and  Plass  (1950)  have  recently  presented  the  view  that  the  major  contribution  to  radiational  cooling 
in  the  atmosphere  is  due  to  the  phenomenon  of  pressure  broadening.  If  one  considers  completely  separated 
absorption  lines,  a certain  amount  of  radiation  from  the  line  wings  will  be  continually  lost  to  space  since  the 
pressure  and  therefore  the  half  width  of  these  lines  decreases  with  height.  This  wing  effect  will  increase  the 
upward  radiation  flux  as  higher  altitudes  are  reached  and  more  and  more  energy,  originating  from  lines  at 
lower  elevations,  penetrates  through  the  atmosphere  without  being  absorbed.  It  is  not  very  likely  that  this 
effect  would  be  of  paramount  importance  in  the  troposphere,  since,  as  was  indicated  in  Section  2,  the  infrared 
absorption  bands  of  both  water  vapor  and  carbon  dioxide  are  composed  of  numerous,  partially  overlapping 
lines.  It  was  indicated  that  for  overlapping  absorption  lines,  wing  radiation  is  probably  reduced  by  square 
root  rather  than  by  linear-pressure  dependence.  Thus  the  major  contribution  to  the  flux  divergence  in  the 
troposphere  must  come  from  the  rapid  decrease  of  the  absorbents  themselves.  In  the  stratosphere,  where 
most  of  the  absorption  lines  become  well  separated,  the  linear-pressure  correction  should  become  valid  and 
the  effect  of  pressure  broadening  could  partially  explain  the  infrared  cooling  above  the  tropopause.  Further 
study  of  the  effect  of  pressure  broadening  on  line  width  and  absorption  is  needed  to  solve  this  problem. 

To  calculate  the  rate  of  radiational  cooling  one  must  be  able  to  determine  the  vertical  divergence  of  the 
net  flux.  If  the  radiation  flux  is  known  at  all  levels,  the  infrared  radiational  temperature  change  per  unit 
time,  in  a unit  volume,  can  be  computed  from  the  first  law  of  thermodynamics, 

Tdd 

dFn  = —cp  — — pdz,  (16) 


where  c,  is  the  specific  heat  of  air  at  constant  pressure  and  8 represents  potential  temperature.  If  we  neglect 
the  effect  of  advection  and  compression,  the  right  side  of  Eq.  (16)  defines  an  effective  temperature  change 
given  by 

dT  i 

dFn  = —cr  — p dz. 
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day-1)  in  equatorial  regions  and  at  about  3 km  (1.0°C  day-1)  at  60°  to  70°  N.  Thus,  even  in  the  case  of  clear 
skies,  the  largest  infrared  cooling  occurs  in  the  middle  troposphere  rather  than  near  the  surface  of  the  earth. 
This  is  essentially  due  to  the  logarithmic  character  of  the  decrease  of  water  vapor  with  height,  particularly 
above  the  turbulence  layer. 

As  a result  of  the  high  surface  temperature  and  very  low  relative  humidity  aloft  in  the  subtropical  region, 
the  layer  near  the  earth’s  surface  loses  heat  more  rapidly  in  the  20°  to  30°  latitude  belt  than  at  any  other 
latitude,  with  the  heat  loss  in  the  surface  layer  being  a minimum  in  the  polar  regions.  Because  of  the  high 
temperatures  in  the  equatorial  regions,  and  the  consequent  high  specific  humidities,  the  largest  rate  of  infra- 
red cooling  aloft,  level  for  level,  is  at  latitudes  0°  to  20°  N.  Here,  the  maximum  of  2.0nC  day-1  is  the  highest 
average  radiative  heat  loss  in  the  atmosphere. 

The  calculations  for  heights  near  the  tropopause  can  only  be  approximate  because  of  our  inadequate 
knowledge  of  the  water  vapor  content  at  these  levels.  The  computed  results  give  an  indicated  rate  of  cooling 
of  0.5°  to  1.0°C  day-1  near  the  tropopause.  These  values  are  in  good  agreement  with  the  results  of  the  calcu- 
lations of  Junge  (1937)  of  the  magnitude  of  the  infrared  cooling  due  to  water  vapor  in  the  lower  stratosphere. 

4.22.  Modifications  Due  to  Average  Cloudiness 

The  cooling  distribution  discussed  above  is  changed  somewhat  if  the  average  cloudiness  for  each  latitude 
belt  is  included  in  the  calculations  for  infrared  radiational  cooling.  The  overall  cooling  of  the  troposphere 
remains  at  about  1.6°C  day-1  in  the  equatorial  region  (0°  to  20°  N).  In  the  polar  region  (50°  to  70°  N), 
however,  as  a result  of  the  predominance  of  low  over  middle  or  high  cloudiness,  the  net  upward  flux  in  the 
lower  layers  is  reduced,  and  the  overall  cooling  rate  is  increased  from  about  0.75°C  day"1  for  clear  skie6  to 
about  1.0°C  day-1  for  conditions  of  average  cloudiness. 

It  is  shown  in  Figs.  11  and  12  that  at  all  latitudes  the  rate  of  infrared  cooling  increases  with  height  from 
the  earth’s  surface  to  about  3 to  5 km  where  there  is  a pronounced  maximum  and  then  decreases  to  the 
tropopause.  By  comparing  the  diagrams  for  clear  and  cloudy  skies,  one  sees  that  clouds  merely  act  as  more 
effective  shields  for  outgoing  radiation  than  does  water  vapor.  Cloudiness  does  not  qualitatively  change  the 
distribution  of  infrared  cooling. 

At  all  latitudes  near  the  earth’s  surface,  the  rate  of  cooling  has  been  decreased  as  a result  of  the  presence 
of  low  clouds.  The  decrease  in  cooling  in  the  lowest  layer  is  about  50  percent  except  in  the  subtropics,  where 
the  amount  of  low  clouds  is  at  a minimum.  Here  the  cooling  has  been  decreased  by  about  30  percent.  An 
additional  effect  of  cloudiness  on  the  cooling  distribution  is  to  define  more  sharply  the  level  of  maximum  cool- 
ing in  the  middle  troposphere.  This  level  of  maximum  cooling  varies  with  the  average  height  of  the  middle 
clouds  and  thus  is  at  a higher  elevation  in  the  tropical  than  in  the  polar  regions.  The  maximum  of  infrared 
cooling  in  the  troposphere  has  been  increased  to  about  2.5°C  day-1  and  is  found  at  latitudes  10°  to  30°  N at 
about  4 to  5 km.  The  average  cloudiness  distribution  does  not  influence  the  infrared  flux  divergence  in  the 
upper  troposphere.  This  is  indicated  by  the  fact  that  the  rate  of  cooling  near  the  tropopause  remains 
essentially  the  same  (0.5°  to  1.0°C  day-1). 

The  heat  lost  by  the  atmosphere  through  long-wave  radiation  is  made  up  in  part  by  the  absorption  of 
insolation  by  water  vapor  and  clouds.  This  absorption  is  relatively  small,  however,  as  will  be  shown  in 
Section  5,  and  the  patterns  of  infrared  cooling  presented  above  dominate  the  distribution  of  net  radiational 
temperature  changes  in  the  atmosphere. 
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4.3.  INFRARED  COOLING  DI  E TO  OZONE  AND  CARBON  DIOXIDE 

In  addition  to  water  vapor,  ozone  and  carbon  dioxide  also  have  absorption  bands  in  the  infrared  region. 
Of  the  three  know  n infrared  absorption  bands  of  ozone,  only  the  band  centered  at  9.6  n might  be  of  meteoro- 
logical significance,  since  it  is  the  strongest  of  the  three.  At  terrestrial  temperatures,  the  maximum  black- 
body  radiation  is  centered  at  about  10  to  12  n,  and  it  is  precisely  in  this  region  in  the  infrared  that  the  water 
vapor  absorption  is  weakest.  One  might  suspect,  therefore,  that  there  is  a significant  amount  of  cooling  due 
to  the  divergence  of  infrared  flux  due  to  ozone.  This  is  found  to  be  the  case  for  the  stratosphere  (30  km). 
(See,  e.g.,  Craig  (1949).)  In  the  troposphere,  however,  the  amount  of  ozone  present  is  extremely  small,  as 
shown  above,  and  both  the  infrared  flux  and  the  divergence  of  the  infrared  flux  due  to  ozone  is  negligible 
compared  with  that  of  water  vapor. 

Carbon  dioxide  has  two  intense  absorption  bands  in  the  infrared,  one  between  12.5  n and  17.5  n and  the 
other  centered  at  about  4.3  /j.  The  former  is  the  more  important  for  atmospheric  radiation,  because  it  is 
both  wider  and  closer  to  the  maximum  black -body  radiation  at  terrestrial  temperatures.  It  has  been  shown 
bv  Callendar  (1941 ) that  for  the  amounts  of  carbon  dioxide  present  in  the  troposphere,  the  net  flux  of  infrared 
radiation  in  the  12.5-  to  17.5-/1  band  is  very  small.  Bruinenberg  (1946)  and  Brooks  (1948)  have  calculated 
the  cooling  due  to  the  flux  divergence  in  the  carbon  dioxide  bands  and  have  found  that  this  cooling  in  the 
lower  and  middle  troposphere  is  of  the  order  of  one  tenth  that  of  water  vapor.  It  has  also  been  suggested  by 
Deacon  (1950)  that  for  layers  thicker  than  100  meters  the  cooling  due  to  carbon  dioxide  is  much  smaller 
than  that  due  to  water  vapor.  For  very  thin  layers  near  the  ground,  however,  the  cooling  could  be  as  large. 

By  the  use  of  the  emissivily  path  length  relationship  for  carbon  dioxide  given  by  Elsasser  (1942),  it  is 
possible  to  compute  directly  the  net  flux  due  to  carbon  dioxide  radiation  at  each  level  in  the  troposphere. 
This  was  done  for  the  tempera  lure-pressure  distribution  of  Table  4 and  the  carbon  dioxide  distribution  of 
Table  5.  The  net  flux  in  the  lower  layers  is  positive  but  very  small,  and  increases  very  slowly  with  height. 
Through  the  troposphere,  the  average  cooling  due  to  carbon  dioxide  is  of  the  order  of  0.1°C  day-1. 

In  the  upper  troposphere  and  in  the  stratosphere,  however,  the  total  amount  of  carbon  dioxide  is  small 
(see  Table  5)  and  the  pressure  correction  to  the  optical  path  for  carbon  dioxide  becomes  important.  The 
cooling  in  the  stratosphere  due  to  carbon  dioxide  is  even  larger  than  that  due  to  ozone  as  has  been  shown  by 
Oder  (1948),  Craig  (1949)  and  Strong  and  Plass  (1950). 

The  infrared  cooling,  due  to  ozone  and  carbon  dioxide,  has  been  omitted  in  the  calculations  of  net 
radiation  temperature  changes  in  the  troposphere.  In  view  of  the  foregoing  discussions,  it  i6  seen  that  this 
omission  affects  only  those  calculations  for  levels  near  the  tropopause. 


5.  CALCULATION  OF  THE  ABSORPTION  OF  INSOLATION  AND  THE  RESULTING  HEATING 

IN  THE  TROPOSPHERE 


5.1.  THE  DISTRIBUTION  OF  INSOLATION  AT  THE  TOP  OF  THE  ATMOSPHERE 

The  incoming  solar  radiation  received  by  a horizontal  surface  of  unit  area  per  unit  time  at  the  top  of  the 
atmosphere  is  given  by 

Io  - (// r3)  cos  f,  (20) 
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Table  9.  The  distribution  of  average  insolation  at  the  top  of  the  atmosphere  during  March. 

(Unit:  cal  cm"*  min  "‘) 


* 

0°  to  10°  N 

10°  to  20°  N 

20“  to  30°  N 

30°  to  40“  N 

40“  to  50“  N 

50“  to  60°  N 

60“  to  70“  N 

I, 

0 6169 

0 5928 

0 5508 

0.4923 

0 4189 

0 3330 

0 2376 

where  f is  the  zenith  angle  of  the  sun,  r is  the  distance  from  the  earth  to  the  sun  in  astronomical  units  (r  — 1.0 
when  the  earth  is  at  mean  distance  from  the  sun)  and  J is  the  average  solar  constant  determined  for  r = 1.0. 
The  value  of  the  average  solar  constant  is  1.94  cal  cm"!  min-1  according  to  Abbot  (1942). 


It  can  be  shown  from  spherical  trigonometric  relationships  (Humphreys  (1940))  that  the  average  inso 
lation  during  any  day  is  given  by 

J 


I o = — (sin  if>  sin  S'wo  + cos  4>  cos  S sin  o>o), 

rr 


(21) 


where  6 is  the  latitude,  5 is  the  declination  of  the  sun  and  u>o  is  the  sun’s  hour  angle  at  sunrise  (or  sunset), 
given  by 

cos  too  = —tan  tan  6.  (22) 

The  average  insolation  for  each  day  during  March  was  calculated  for  each  latitude  belt  by  use  of  Eqs. 
(21)  and  (22).  Values  of  r2  and  6 for  each  day  are  given  in  the  American  Ephcmeris.  The  average  of  the 
values  of  r2  and  6 for  the  four  years  1920, 1925,  1930  and  1935  was  used  to  eliminate  the  inequalities  of  5 and  r 
on  the  same  date  in  different  years  caused  by  the  fact  that  a year  does  not  contain  an  integer  number  of  days. 

In  calculating  the  average  values  of  h for  March,  the  arithmetic  mean  was  taken  of  the  daily  average 
values  during  the  month.  The  average  values  of  l0  were  computed  at  latitudes  5°,  15°  . . . 65°  to  approxi- 
mate the  means  for  the  10°  latitude  belts  0°  to  70°  N and  are  shown  in  Table  9. 


5.2.  THE  TOTAL  ABSORPTION  OF  SOLAR  ENERGY  BY  MOLECULAR  OXYGEN 

The  solar  beam  which  enters  the  top  of  the  atmosphere  is  only  slightly  depleted  before  reaching  the 
ozonosphere.  This  depletion  is  to  a large  extent  due  to  the  absorption  characteristics  of  molecular  oxygen. 

Penndorf  (1950)  has  given  an  estimate  of  the  vertical  distribution  of  the  energy  absorbed  per  unit 
volume  per  unit  time  under  the  assumptions  that: 

1.  absorption  by  molecular  oAygen  takes  place  in  the  wavelength  region  0.11  to  0.22  ju, 

2.  insolation  in  this  wavelength  interval  is  at  an  effective  black-body  temperature  of  5300°K, 

3.  zenith  distance  is  0°,  and 

4.  the  rate  constant  of  the  3-body  collision  producing  molecular  oxygen  in  the  E-layer  is 
ki  = 1.3  X 10“ **  t/To  cm*  sec-1. 

The  total  energy  absorbed  per  unit  volume  per  unit  time  is  given  by 

dQ/dv  dt  = he  I Ax  X-1  d\ 

J*  i 

where  Xi  and  X2  are  the  wavelength  limits  of  the  integration,  A\  is  the  number  of  quanta  absorbed  per  cm1, 
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h is  Planck’s  constant  and  c is  the  speed  of  light.  Penndorrs  results  show  absorption  of  10~u  cal  cm-3  sec-1 
at  about  90  km  increasing  to  a maximum  at  102  km  of  about  1.5  X 10“ 11  cal  cm"'  sec-1  and  then  decreasing 
to  about  10_li  cal  cm-'  sec-1  at  about  125  km. 

Results  of  spectrograph  measurements  made  during  rocket  ascents  have  shown  that  at  heights  of  50  km, 
there  is  still  complete  absorption  below  about  0.23  m*  If  we  assume  constant  absorption  from  50  to  90  km, 
the  total  absorption  in  a unit  column  of  air  per  unit  time  would  be 

/•14* 

dQ/dt  ~ / (dQ/dv  dt)  ■ dz 

J 60 

^ 10“ 4 cal  cm“!  min-1. 

The  assumptions  made  above,  upon  which  this  value  is  based,  would  give  the  probable  maximum  absorp- 
tion by  molecular  oxygen.  This  value  is  smaller  than  the  total  ozone  absorption  by  a factor  of  200. 


5.3.  THE  ABSORPTION  OF  INSOLATION  BY  ATMOSPHERIC  OZONE 

It  has  already  been  indicated  that  the  direct  solar  beam  remains  almost  completely  undiminished  until 
it  reaches  the  top  of  the  ozone  layer.  Through  the  ozonosphere,  however,  most  of  the  energy  in  the  ultra- 
violet region  of  the  solar  spectrum,  as  well  as  in  a small  portion  of  the  visible  part  of  the  spectrum,  is  absorbed. 
The  energy  absorbed  by  ozone  in  the  atmosphere  depends  on  the  distribution  of  energy  in  the  6olar  spectrum, 
the  absorption  coefficients  of  ozone,  and  the  total  amount  and  vertical  distribution  of  the  ozone  in  a column 
traversed  by  the  solar  beam. 

The  energy  distribution  of  the  solar  spectrum  is  very  nearly  that  of  a black-body  radiator  at  an  effective 
temperature  of  about  6000°K.  Measurements  by  Abbot  et  al.  (1923)  and  Pettit  (1932,  1940),  among 
others,  gave  indications  that  below  0.4  m the  solar  energy  curve  drops  sharply  from  that  corresponding  to  a 
black  body  at  6000°K.  These  results  have  been  confirmed  by  more  recent  observations  taken  during  rocket 
flights  at  White  Sands,  New  Mexico.  The  rocket  measurements  indicate  that  the  insolation  at  the  top  of 
the  atmosphere  ranges  from  about  80  percent  at  0.4  m to  about  40  percent  at  0.3  m,  and  to  about  10  to  15  per- 
cent at  0.25  n,  of  that  of  a black -body  radiator  at  6000°K  (e.g.,  Hulbert  (1947),  Newell  (1950)). 

Ozone  has  absorption  bands  in  the  ultraviolet,  visible  and  infrared  regions  of  the  spectrum.  The  most 
important  bands  for  the  depletion  of  insolation  are  the  Hartley  (0.22  to  0.32  n),  the  Huggins  (0.32  to  0.36  n) 
and  the  Chappuis  (0.44  to  0.76  m)  bands.  The  absorption  coefficients  in  these  regions  have  been  the  subject 
of  many  extensive  investigations,  the  most  important  of  which  were  those  conducted  by  Colange  (1927), 
Ny  and  Choong  (1932,  1933)  and  Vassey  (1941)  whose  results  seem  to  be  the  most  consistent  and  reliable. 
Craig  (1950)  has  recently  presented  an  excellent  summary  of  the  extant  data  concerning  the  absorption 
coefficients  of  ozone  in  these  three  bands.  The  data  show  very  intense  absorption  in  the  Hartley  band,  with 
an  absorption  coefficient  at  the  band  center  (about  0.25  m)  of  about  6 X 102  cm-1  03  at  NTP,  decreasing 
absorption  in  the  Huggins  band  with  an  average  absorption  coefficient  of  about  10“'  cm-1  03  at  NTP,  then 
increasing  again  through  the  Chappuis  band  to  a maximum  a t about  0.59  m of  slightly  larger  than  10-1  cm-1  Oj 
at  NTP.  The  absorption  bands  for  ozone  in  the  infrared  region  are  centered  at  4.75  n,  9.6  m and  14.1  m.  They 
contribute  very  little  to  the  depletion  of  incoming  radiation,  since  the  insolational  energy  in  the  infrared 
region  is  very  small. 
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The  monochromatic  depletion  of  insolation  by  ozone  in  a horizontal  layer  of  vertical  depth  dz  is  given  by 

d/x  - I\kx  n m dz,  (23) 

where  I\  is  the  monochromatic  intensity  of  a parallel  beam  at  the  top  of  the  layer,  k\  is  the  absorption  coeffi- 
cient of  ozone  at  wavelength  X,  n is  the  ozone  density  at  height  z and  m is  the  relative  length  of  the  path 
through  which  the  solar  beam  passes  and  depends  on  the  zenith  angle.  If  Eq.  (23)  is  integrated  from  height  z 
to  the  top  of  the  atmosphere,  we  get 


h = /.>  exp(-  k\n  m dz). 


(24) 


Upon  substitution  of  Eq.  (24)  into  Eq.  (23)  we  get  an  expression  for  the  energy  absorbed  per  unit  volume 
along  the  solar  beam.  The  energy  absorbed  in  a wavelength  interval  dX  is  then 

dh 


m dz 


dX  = Io,xk\it  exp(  — fcx  mN)d\, 


(25) 


where 


N 


■/: 


n(z)  dz. 


The  distribution  of  n as  a function  of  height  and  latitude  during  March  has  been  given  in  Section  2. 


Since  the  integration  of  Eq.  (25)  need  be  carried  out  only  over  the  regions  of  significant  ozone  absorption 
(i.e.,  where  /ox  kx  is  not  negligible)  it  is  important  to  consider  only  the  absorption  in  the  Hartley,  Huggins  and 
Chappuis  bands.  For  fixed  values  of  N and  m,  Eq.  (25)  can  be  integrated  numerically.  Craig  (1949)  has 
presented  the  results  of  such  an  integration  in  the  form  of  a graph  of  absorption  per  unit  volume  per  unit 
ozone  density  against  total  path  length  of  ozone.  For  the  relative  path  length  m,  Craig  used  sec  f.  This  is 
certainly  valid  for  f < 75°.  For  f > 75®  the  absorption  is  small  since  /#  is  small,  and  at  f « 90°  the  absorp- 
tion is  zero  because  J*  = 0. 


The  total  optical  path  and,  therefore,  the  absorption  by  ozone,  changes  considerably  during  the  day. 
The  average  absorption  is  given  by 


a. 


(26) 


where  w is  the  hour  angle  (o>  » o>0  at  sunrise  or  sunset),  a,„  is  the  absorption  per  unit  volume  at  hour  angle  u 
and  height  z.  For  constant  u>  there  is  little  variation  in  N sec  f from  day  to  day.  Thus,  values  of  N sec  f 
for  March  16  were  used  to  represent  the  average  ozone  path  length  during  March.  Values  of  absorption  per 
unit  volume  per  unit  ozone  density  for  each  hour  from  sunrise  to  noon  were  obtained  as  a function  of  iV  sec  / 
from  the  chart  prepared  by  Craig  (1949).  The  integral  in  Eq.  (26)  was  then  approximated  by  use  of 
Weddle’s  rule  where  <■>«  ft;  w/2.  Thus  the  average  absorption  during  the  day  is. 


a«,  1/ 40[otio  + 5otii  + o.i  -f  6of>»  + o,,«  + 5o,,»  + o,t«J 


(27) 


where  the  subscript  a,,,-  (i  * 1 ...  6)  refers  to  the  absorption  per  unit  volume  at  height  z and 
hour  angle  iir/12. 

The  distribution  of  absorption  per  unit  volume  by  atmospheric  ozone  is  shown  in  Fig.  13.  The  absorp- 
tion is  rather  large  in  the  upper  ozonosphere  (40  to  50  km)  where  the  principal  absorption  is  in  the  Hartley 
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Fig.  13.  Distribution  of  ab- 
sorption of  insolation  (10** 
cal  cm*2  min*1  km*1)  by  at- 
mospheric ozone. 
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band  and  increases  to  a maximum  of  about  5 X 10-4  cal  cm-3  min-1  km-1  at  about  25  to  30  km  where  the 
remaining  insolation  in  the  Hartley  band  is  completely  absorbed.  Below  about  15  to  20  km  (through  the 
troposphere),  the  absorption  is  relatively  small,- being  about  0.75  X 10-4  cal  cm-1  min-1  km-1,  and  is  due 
almost  entirely  to  absorption  in  the  weaker  Huggins  and  Chappuis  bands.  The  total  absorption  per  unit 
area  increases  rapidly  through  the  ozonosphere  and  is  practically  constant  in  the  troposphere  at  a value  of 
about  0.01  cal  cm'1  min-1.  This  total  absorption  is  almost  unchanged  with  latitude  from  0°  to  70°  N, 
resulting  from  the  fact  that  the  insolation  on  a unit  horizontal  surface  decreases  while  the  total  ozone  content 
increases  with  increasing  latitude.  North  of  70°,  it  appears  that  the  total  absorption  should  decrease 
rapidly,  since  both  the  total  ozone  content  and  the  vertical  component  of  insolation  decrease. 
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The  value  of  0.01  cal  cm"2  min  1 represent*  about  1.6  percent  of  the  incoming  radiation  at  0°  to  10°N, 
increasing  to  about  4 percent  of  the  insolation  at  60s  to  70°  IS.  The  average  absorption  in  the  atmosphere 
due  to  ozone  is  about  2.1  percent  of  the  average  insolation  for  the  northern  hemisphere  during  March.  This 
value  of  percentual  absorption  is  lower  than  that  given  by  Dobson  (1942)  of  5 to  8 percent  and  equal  to  that 
given  by  Fritz  (1949)  of  2.1  percent  for  the  mean  absorption  by  ozone  for  the  year.  The  higher  value  is 
based  on  the  assumption  of  a 6000°K  black -body  solar  energy  curve.  It  has  already  been  pointed  out,  how- 
ever, that  the  ultraviolet  insolation  received  at  the  top  of  the  atmosphere  is  only  about  40  percent  of  that  of  a 
black -body  radiator  at  6000°K. 

The  absorption  per  unit  volume  in  the  troposphere  by  atmospheric  ozone  is  exceedingly  small.  The 
contribution  to  heating  in  the  troposphere  due  to  this  absorption  is  about  0.02eC  day-1  near  the  tropopause 
decreasing  to  about  0.004°C  day-1  near  the  earth’s  surface.  This  represents  a minor  component  of  heating 
in  the  troposphere  and  is  neglected  in  the  present  calculations  of  net  radiationai  temperature  changes. 


5.4.  THE  ABSORPTION  BY  WATER  VAPOR 

The  absorption  of  insolation  by  water  vapor  is  almost  completely  confined  to  the  visible  and  near 
infrared  regions  since  about  99  percent  of  the  energy  of  the  solar  spectrum  is  contained  in  wavelengths  below 
4/i.  This  absorption  has  been  studied  experimentally  by  Fowle  (1915,  1917)  and  llettner  (1918),  among 
others.  Fowle  determined  the  fractional  absorption  of  light  from  a Nernst  glower  as  a function  of  water 
vapor  path  length.  The  light  was  sent  through  a long  tube  containing  known  concentrations  of  water  vapor 
at  room  temperature  and  atmospheric  pressure  and  the  absorption  in  each  water  vapor  band  was  determined 
spectroscopically.  The  total  absorption  of  insolation  as  a function  of  total  water  vapor  was  then  calculated 
by  correcting  for  the  difference  between  the  strength  of  the  light  source  and  the  average  solar  constant  of 
1.94  cal  era-2  min-1.  Hettner’s  measurements,  on  the  other  hand,  were  carried  out  for  water  vapor  at  a 
temperature  of  400°K  and  are  not  directly  applicable  to  meteorological  studies.  The  results  given  by 
Hettner  for  absorption  values  in  the  near  infrared  are  much  larger  than  are  found  for  atmospheric  water 
vapor  (see  Yamamoto  and  Onishi  (1949  a)). 

The  absorption  of  insolation  by  water  vapor  in  the  atmosphere  has  been  determined  directly  by  Kimball 
(1927,  1930)  and  Hoelper  (1937,  1943).  These  studies,  based  on  solar  radiation  measurements,  give  the 
distribution  of  absorption  as  a function  of  total  precipitable  water  in  the  path  of  the  solar  beam.  The  dis- 
tribution of  water  vapor  was  found  from  various  empirical  formulas.  A table,  summarizing  the  results  of 
Fowle,  Kimball,  and  Hoelper,  is  given  by  Hoelper  (1943)  where  a comparison  of  the  various  results  shows  that 
Fowle’s  corrected  values  for  absorption  of  insolation  are,  for  all  water  vapor  path  lengths,  about  10  to  25 
percent  higher  than  the  values  given  by  Kimball  and  Hoelper,  the  latter  two  sets  of  values  being  in  good 
agreement  with  each  other. 

Formulas  for  the  absorption  per  unit  volume  at  any  height  with  any  zenith  distance  of  the  sun  were 
obtained  by  Yamamoto  and  Onishi  (1951)  using  the  theoretical  values  of  the  water  vapor  absorption  coeffi- 
cients obtained  by  them.  Their  results  (which  include  a pressure  correction)  agree  with  those  of  Kimball 
and  Hoelper. 

Inherent  in  the  observations  used  by  Kimball  and  Hoelper  was  the  effect  of  pressure  broadening  on  the 
near  infrared  water  vapor  absorption  bands.  If  the  total  precipitable  water  in  the  atmosphere  is  replaced 


u 


by  the  effective  water  vapor  path  length,  taking  into  account  a square-root  pressure  correction.  Fowl#'* 
results  are  reduced  to  the  average  of  the  values  given  by  Kimball  and  Hoelper. 

Miigge  and  M oiler  (1932)  have  shown  that  Fowle’s  values  can  be  fitted,  to  a fair  degree  of  accuracy,  by 
the  empirical  formula 

a = 0.172  («)••*,  (28) 

where  a represents  the  energy  absorbed  in  cal  cm-1  min-1  in  a column  parallel  to  the  solar  beam  containing 
u gm  cm'*  of  precipitable  water.  If  u is  measured  in  a vertical  column,  and  we  replace  u by  u*  (the  effective 
optical  depth  for  water  vapor)  the  absorption  along  the  solar  beam  becomes 

a = 0.172(u*  sec  f)nJ,  (29) 

where,  as  in  previous  discussions, 

u*  *■  Vp/po  u. 

The  average  of  the  values  given  by  Kimball  and  by  Hoelper  can  be  expressed  by  a relationship  analogous 
to  Eq.  (29).  That  is, 

a = 0.147(u  sec  f)"-»,  (30) 

the  difference  being  that  no  pressure  correction  is  used  in  Eq.  (30)  for  the  water  vapor  distribution.  The  two 
expressions  for  the  water  vapor  absorption  given  above  give  almost  the  same  values  for  a wide  range  of  water 
vapor  path  lengths.  For  convenience,  since  the  square-root  pressure  correction  was  used  in  the  determination 
of  total  precipitable  water  as  given  in  Section  3,  the  approximation  given  in  Eq.  (29)  was  used  in  the  present 
calculations. 

The  depletion  of  insolation  by  water  vapor  in  a vertical  column  extending  from  height  z to  the  top  of  the 
atmosphere  is  then  given  by 

a,j  — 0. 172(u*(z)  sec  f]0J  cos  f.  (31) 

For  the  average  absorption  during  the  day,  we  can  write  an  analogous  expression  to  that  used  in  the  case  of 
absorption  by  ozone.  That  is, 

a,  = l/2«o  a,f  do>,  (32) 

where  the  integral  in  Eq.  (32)  is  again  computed  numerically  by  the  use  of  Weddle’s  rule. 

As  was  the  case  for  ozone  for  constant  u,  a,  [ varies  very  little  from  day  to  day.  To  avoid  a great 
number  of  tedious  computations,  a,  as  calculated  for  March  16  was  used  to  represent  average  values  during 
March.  Values  of  5,  were  computed  from  Eqs.  (31)  and  (32)  for  each  10*  latitude  belt  (0°  to  70°  N)  for 
each  kilometer  of  height  in  the  troposphere.  The  vertical  distribution  of  absorption  during  March  for  each 
latitude  belt  is  given  in  Table  10.  The  absorption  by  water  vapor  in  a unit  vertical  column  extending  from 
the  surface  to  the  top  of  the  atmosphere  amounts  to  about  0.093  cal  cm-1  min-1  at  0°  to  10°  N decreasing  to 
about  0.021  cal  cm'*  min-1  at  60°  to  70°  N.  This  represents  an  absorption  of  about  IS  percent  of  the  insola- 
tional  energy  received  at  the  top  of  the  atmosphere  in  equatorial  regions  decreasing  to  about  9 percent  of  that 
received  in  the  region  60°  to  70°  N. 

The  distribution  of  absorption  is  shown  by  a series  of  curves  in  Fig.  14.  These  curves  are  drawn  for 
each  10°  latitude  belt  from  the  surface  to  10  km.  The  variation  of  absorption  with  height  can  be  represented 
with  good  accuracy  by  a relationship  of  the  form  a = o0  exp  {—bz)  where  the  coefficient  b varies  from  about 
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Fig.  14.  Vertical  variation  of 
absorption  of  insolation  by 
water  vapor  for  each  10°  lati- 
ta4e  bolt,  cUar  tkUt. 


Table  10.  The  average  absorption  of  mentation  by  atmospheric  water  vapor  daring  March. 

(Unit:  eel  ear*  mh>->) 


*(km) 

0*  to  10*  N 

10°  to  20°  N 

20°  to  30°  N 

30*  to  40°  N 

40°  to  50°  N 

50°  to  60°  N 

60*  to  70*  N 

0 

0.0929 

0.0839 

0.0735 

0.0577 

0.0458 

0.0324 

0.0206 

1 

0.0796 

0.0713 

0.0608 

0.0496 

0.0402 

0.0279 

0.0182 

a 

0.0684 

0.0596 

0.0500 

0.0421 

0.0344 

0.0238 

0.0157 

3 

O.OS86 

0.0487 

0.0416 

0.0344 

0.0279 

0.0198 

0.0133 

4 

0.0498 

0.0398 

0.0346 

0.0289 

0.0224 

0.0162 

0.0112 

5 

0.0416 

0.0329 

0.0283 

0.0242 

0.0182 

0.0134 

0.0098 

6 

0.0344 

0.0267 

0.0229 

0.0193 

0.0148 

0.0113 

0.0087 

7 

0.0281 

0.0213 

0.0183 

0.0152 

0.0122 

0.0098 

0.0079 

8 

0.0227 

0.0169 

0.0146 

0.0124 

0.0105 

0.0087 

0.0071 

9 

0.0180 

0.0134 

0.0116 

0.0106 

0.0095 

0.0079 

0.0065 

10 

0.0141 

0.0107 

0.0094 

0.0094 

0.0088 

0.0072 

0.0060 

11 

0.0108 

0.0086 

0.0079 

0.0086 

0.0082 

12 

0.0083 

0.0070 

0.0067 

0.0080 

13 

0.0064 

0.0060 

0.0060 

0.0074 

14 

0 0055 

0.0053 

0.0056 

IS 

0.0048 

0.0048 

0.0051 

16 

0 0044 

0 0044 

0.0047 

0.16  at  low  latitudes  to  about  0.13  in  the  more  northerly  regions.  This  form  of  the  absorption  curves  follows 
from  the  essentially  logarithmic  character  of  the  effective  water  vapor  distribution  (see  Section  3). 


5.5.  HEATING  DUE  TO  WATER  VAPOR  ABSORPTION 

The  average  heating  in  a thin  layer  due  to  the  absorption  of  insolation  by  water  vapor  can  be  found  by 
computing  the  absorption  in  the  layer  du.  Tanck  (1940)  has  given  a formula,  baaed  on  the  empirical  expres- 
sion of  Mugge  and  Moller,  for  the  direct  computation  of  average  heating  daring  a solar  day  resulting  from 
absorption  of  insolation  by  atmospheric  water  vapor.  The  average  temperature  change  in  units  of  °C  min*1 
was  given  as 

dT  .0521  q 
it  “ (u)Mc, 


(33) 
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Fig.  15.  Distribution  of  beating 
in  the  troposphere  (°C  day*1) 
doe  to  the  absorption  of  in- 
solation by  water  vapor,  clear 
skies. 


where  q is  the  average  specific  humidity  in  the  layer,  and  Z is  the  average  value  of  cos  °-7f  integrated  over 
the  solar  day. 

If  the  vertical  distribution  of  absorption  is  known,  however,  the  temperature  changes  in  finite  layers  due 
to  absorption  of  insolation  can  be  computed  with  sufficient  accuracy  by  the  method  used  to  determine  infra- 
red cooling.  Thus, 

AT  * 5.9  X 10*  Aa/Ap,  (34) 

where  AT  is  the  temperature  change  in  °C  day-1  due  to  the  absorption  in  a layer  of  A p mb  thick,  if  compres- 
sions! effects  are  neglected. 

5.51.  Clear  Skies 

Temperature  changes  resulting  from  the  absorption  of  insolation  in  the  troposphere  have  been  investi- 
gated by  Moller  (1935),  Tanck  (1940)  and  others.  By  the  use  of  Mugge  and  M oiler’s  adaptation  of  Fowle’s 
absorption  distribution  (discussed  above),  Moller  determined  the  average  vertical  distribution  of  beating 
for  June  over  Lindenberg  (52°  N).  His  results  show  a minimum  of  heating  at  1 km  of  about  0.2°C  day-1 
increasing  to  a maximum  of  about  0.5°C  day-1  at  6 km  and  then  decreasing  above  6 km  Tanck,  using 
Eq.  (33),  calculated  the  effect  of  the  temperature  and  relative  humidity  distribution  on  heating  in  the 
atmosphere  by  water  vapor  absorption.  Two  different  temperature  curves  were  assumed,  and  for  each  of 
these  he  assumed  three  relative  humidity  distributions  corresponding  to  conditions  in  the  center  of  an 
anticyclone  and  directly  east  and  west  of  a cyclone.  In  all  cases  where  the  relative  humidity  increased  or 
remained  constant,  the  maximum  warming  (about  0.7°C  day-1)  occurred  at  the  top  of  the  layer  of  constant 
or  increasing  relative  humidity.  He  also  computed  the  heating  distribution  up  to  6 km  for  data  representing 
a three-year  average  (1936-1938)  of  temperature  and  moisture  over  Hamburg  (53°  N).  The  average  curves 
for  March  show  a rate  of  heating  of  about  0.6°C  day-'  from  the  surface  to  about  3 km  decreasing  to  about 
0.4°C  day-1  at  6 km. 
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Iq  the  present  calculations  the  distribution  of  heating  in  the  troposphere  by  the  absorption  of  water 
vapor  was  determined  for  each  10°  latitude  belt.  The  temperature  changes  were  computed  by  means  of 
Eq.  (34)  and  Table  10  and  therefore  include  the  effect  of  pressure  on  the  distribution  of  absorption.  The 
variation  of  heating  with  latitude  and  with  height  is  shown  in  Fig.  15  for  clear  skies. 

The  results  generally  confirm  the  order  of  magnitude  found  by  previous  investigators.  Although  small 
when  compared  to  the  rate  of  infrared  cooling,  the  values  of  the  computed  heating  are  by  no  means  negligible, 
amounting  in  some  cases  to  about  50  percent  of  the  temperature  changes  due  to  infrared  cooling.  The 
largest  heating  values  (about  0.7°C  day-1)  are  found  in  the  lower  troposphere  (0  to  5 km)  in  the  region 
0°  to  20°  N.  Despite  the  longer  path  lengths  of  the  solar  beam,  the  heating  decreases  northward  (0.1°C 
day-1  at  60°  to  70°  N)  because  of  the  rapid  decrease  of  total  precipitable  water  with  latitude.  In  the  upper 
troposphere  (at  10  km)  the  heating  is  about  0.4°C  day-1  at  0°  to  10°  N decreasing  to  about  0.1°C  day-1  at 
60°  to  70°  N.  This  results  from  the  fact  that  even  at  10  km  our  physical  model  shows  more  precipitable 
water  in  the  equatorial  than  in  polar  regions.  It  should  be  pointed  out  that  the  heating  values  in  the  upper 
troposphere  are,  at  best,  merely  estimates  of  the  magnitude  of  the  actual  temperature  changes. 

5.52.  Reflection,  Transmission  and  Absorption  by  Clouds 

The  distribution  of  heating  discussed  in  the  preceding  section  is  considerably  altered  by  the  presence  of 
clouds  in  the  atmosphere.  This  results  from  the  fact  that  absorption  of  insolation  by  cloud  particles  in  the 
near  infrared  is  somewhat  larger  than  absorption  by  water  vapor  and,  more  important,  the  large  scattering 
by  water  droplets  greatly  reduces  the  transmitted  insolation  that  could  be  absorbed  by  water  vapor  below 
the  cloud. 

The  generally  accepted  value  of  the  average  reflectivity  of  clouds  has  heretofore  been  that  given  by 
Aldrich  (1919).  Aldrich’s  value  of  the  cloud  albedo  (78  percent)  was  based  on  the  average  of  about  100  ob- 
servations of  the  incident  and  reflected  radiation  from  the  tops  of  stratus  and  fog  in  the  valley  adjacent  to 
Mt.  Wilson.  No  apparent  relationship  was  found  between  cloud  thickness  and  reflectivity. 

Neiburger  (1949)  has  recently  presented  observations  of  the  absorption  and  reflection  of  clouds  under 
conditions  similar  to  those  studied  by  Aldrich.  Almost  simultaneous  pyrheliometric  measurements  of  total 
upward  and  downward  radiation  were  made  above  and  below  stratus  clouds,  on  about  50  days  during  the 
summer  of  1945.  He  found,  contrary  to  Aldrich’s  conclusions,  that  the  reflectivity  changed  markedly  with 
changes  in  cloud  depth,  being  very  small  when  the  cloud  is  thin  and  increasing  to  about  80  percent  when  the 
cloud  depth  is  about  2000  ft.  The  average  distribution  of  drop  size  in  the  clouds  studied  shows  a mode  at  a 
radius  of  about  7 X 10~4  cm.  The  liquid  water  content  was  found  to  increase  from  the  base  to  the  top  of  the 
clouds  with  an  average  water  density  of  about  0.3  gm  m~*. 

The  effect  of  the  size  and  density  of  cloud  particles  and  the  thickness  of  the  cloud  on  the  reflection,  trans- 
mission and  absorption  of  solar  radiation  by  clouds  has  been  studied  theoretically  by  Hewson  (1943).  Hew- 
son  showed  that  for  a cloud  with  droplet  radius  of  5 X 10~4  cm  and  liquid  water  density  1.0  gm  m-*,  the 
reflectivity  increases  with  cloud  depth  from  about  38  percent  for  a depth  of  20  meters  to  about  93  percent  for 
a depth  of  4000  meters.  The  absorptivity  by  the  cloud,  under  the  same  conditions  of  drop  size  and  density, 
increases  from  0.6  percent  at  20  meters  to  about  6.5  percent  for  a depth  of  4000  meters.  It  was  also  shown 
that  for  a cloud  200  meters  thick,  80  percent  of  the  absorption  takes  place  in  the  top  100  meters  of  the  cloud. 
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The  results  published  by  Hewson  indicate  that  there  is  practically  no  variation  of  reflection,  transmission 
or  absorption  with  solar  zenith  distance.  Thig  is  generally  confirmed  by  the  data  published  by  Haurwitz 
(1948)  of  the  ratio  between  insolation  received  on  a horizontal  surface  at  Blue  Hill  Observatory  for  overcast 
skies  of  various  cloud  types  to  that  received  for  clear  skies.  The  data  indicates  that  the  transmission  decreases 
from  about  80  percent  for  high  clouds  (Ci  to  Cs)  to  about  20  percent  for  fog  and  Ns.  There  is  very  little 
change  in  transmission  shown  for  a change  in  zenith  angle  from  0°  to  75°.  Since  all  of  the  cloud  observations 
were  made  at  the  ground,  there  were  no  reports  of  cloud  thickness  to  relate  the  change  of  transmission  to 
cloud  depths. 

From  the  results  discussed  above,  it  is  evident  that  Aldrich’s  value  for  the  albedo  of  clouds  is  probably 
correct  for  typical  stratus  about  2000  ft.  thick,  but  does  not  represent  the  average  reflectivity  of  clouds  inde- 
pendent of  cloud  structure  and  cloud  depth.  The  distribution  of  percentual  reflection  and  transmission  of 
incident  insolation  as  a function  of  cloud  depth,  as  determined  by  Hewson,  is  shown  in  Fig.  16  for  a drop 
radius  of  5 X 10-4  cm.  Three  sets  of  curves  are  drawn  corresponding  to  liquid  water  densities  of  0.1,  1.0 
and  5.0  gm  m~*.  The  absorption  distribution  for  a constant  drop  size  and  cloud  density  is  determined  from 
the  curves  of  reflection  and  transmission  since  the  percentual  absorption  is 

A - 1 - (R  -|-  P),  (35) 

where  R is  the  reflectivity  and  P is  the  transmissivity.  A set  of  curves  based  on  the  data  taken  from  Nei- 
burger’s  paper  is  included  to  show  the  consistency  between  the  observed  and  theoretically  derived  data. 

It  is  seen  from  the  curves  in  Fig.  16,  that  for  constant  drop  size  and  cloud  density,  the  absorption  and 
reflection  increase  while  the  transmission  decreases  with  increasing  cloud  depth. 

There  is,  of  course,  a wide  variation  in  drop  size  and  cloud  density  not  only  from  cloud  to  cloud,  but  also 
within  the  same  cloud.  For  the  present  investigation  it  is  sufficient  to  consider  average  cloud  structures  by 
cloud  type  and  to  assume  approximate  homogeneity  within  each  cloud  group.  Hewson  (1943)  quotes 
Bricard’s  measured  values  of  drop  radii  for  various  types  of  clouds.  Similar  values  are  given  by  aufm  Kampe 
(1950)  based  on  the  airplane  measurements  of  Diem.  Aufm  Kampe  also  gives  values  of  the  liquid  water 
content  computed  for  various  types  of  clouds. 

Table  11  presents  the  values  of  average  drop  radius  and  liquid  water  content  for  each  of  the  cloud  groups 
used  in  the  present  study. 


Table  11.  Average  drop  radius  and  liquid  water  density  by  cloud  groups. 


Cloud  group 


Radius  (10"'  cm)  Density  gm  m~* 


Cirrus 

Cumulonimbus 

Altoetratus-altocumulus  . . . 

Cumulus 

Stratocumulus-fractocumulus 
Nimboetratus 


1 

10 

5 

5 

5 

10 


0 1 
2 5 
0 2 
05 
0.3 
0 5 


The  cloud  depths,  in  general,  vary  with  latitude.  The  thicknesses  assumed  in  the  following  discussion  for 
each  cloud  group,  are  based  on  the  values  of  the  tops  and  bases  given  for  each  latitude  belt  in  Section  3 of  this 
report.  For  the  As  and  Stcu  cloud  depths,  a constant  thickness  of  200  meters  was  used  for  all  latitudes. 
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Fig.  17.  Distribution  of  heating 
in  the  troposphere  (°C  day*1) 
doe  to  the  absorption  of  in- 
solation by  water  vapor  and 
clonda,  average  cloudiness. 


5.53.  Absorption  and  Heating  for  Average  Cloudiness 

It  has  been  pointed  out  previously  that  the  heating  in  the  atmosphere,  resulting  from  absorption  of  inso- 
lation, is  affected  by  tbe  presence  of  clouds  because  of  the  absorption  of  the  cloud  layer  and  the  reduced  trans- 
mission of  insolation  resulting  from  the  high  reflectivity  of  the  clouds. 

The  inclusion  of  absorption  by  clouds  in  the  calculations  of  heating  in  the  troposphere  introduces  a num- 
ber of  factors  which  increase  the  complexity  of  the  computations  out  of  proportion  to  the  increased  accuracy 
of  the  heating  values.  The  following  assumptions  were  therefore  made  in  order  to  simplify  the  computations 
and  still  maintain  a reasonable  picture  of  tbe  role  played  by  clouds  in  tbe  heating  of  the  atmosphere. 

1.  There  is  no  overlap  of  clouds  along  the  solar  beam. 

2.  The  absorption  by  the  cloud  layer  contributes  to  the  total  absorption  in  the  kilometer  interval  of 
which  the  cloud  top  is  a part. 

3.  Absorption  at  the  sides  of  clouds  is  small  and  can  therefore  be  omitted  in  tbe  calculation. 

4.  Absorption  by  water  vapor  of  the  insolation  reflected  by  clouds  and  tbe  earth’s  surface  is  small 

and  can  be  neglected. 

• 

The  percentage  absorption  by  clouds  for  each  latitude  belt  was  computed  for  the  thickness  of  each  cbud 
group  from  the  appropriate  curves  in  Fig.  16.  The  energy  absorbed  by  the  cloud  layer  in  cal  cm-*  min-1  was 
then  found  by  considering  the  insolation  incident  to  the  cloud  and  the  percentage  of  sky  cover  for  each  cloud 
group.  To  this  was  added  the  energy  absorbed  by  Water  vapor  for  each  kilometer  interval.  The  latter  was 
found  by  reducing  the  absorption  by  water  vapor  determined  in  the  case  of  clear  skies  by  a factor  based  on  the 
reduced  transmission  of  insolation  resulting  from  the  presence  of  clouds. 
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Tbe  distribution  of  beating  in  the  tropoephere  considering  average  cloudiness  is  calculated  from  Eq.  (34) 
where  Ap  represents  the  pressure  difference  for  each  kilometer  interval.  This  distribution  (Fig.  17)  shows 
that  the  m«vin»nn  heating  is  no  longer  found  near  the  surface  but  at  the  level  of  the  middle  clouds.  In  the 
equatorial  regions  the  largest  heating  values  are  at  about  5 km  (0.8°C  day-1)  and  in  the  polar  regions  the 
largeet  heating  values  are  at  about  3.5  km  (0.25°C  day-1).  It  is  interesting  to  note  that  although  the  local 
heating  at  the  predominant  cloud  levels  has  been  increased  as  a result  of  the  absorption  by  clouds,  the  heating 
below  the  clouds  has  been  reduced.  This  is  due  to  the  fact  that  the  clouds  transmit,  on  the  average,  only 
about  60  to  70  percent  of  the  incident  insolation.  The  reduced  transmission  results  in  a noticeable  reduction 
in  the  absorption  by  water  vapor  below  tbe  cloud  layers.  The  heating  of  the  upper  troposphere  is  not  affected 
by  cloudiness.  This  is  borne  out  by  a comparison  of  Figs.  15  and  17.  Just  below  the  tropopause  the  temper- 
ature changes  shown  on  both  charts  are  very  small,  amounting  to  about  0.1  to  0.2°C  day*1. 


6.  THE  DISTRIBUTION  OF  NET  RADIATIONAL  TEMPERATURE  CHANGE 

IN  THE  TROPOSPHERE 

The  net  radiational  temperature  change  in  the  tropoephere  is  the  difference  between  tbe  infrared  cooling 
and  the  heating  due  to  absorption  of  insolation.  For  both  clear  and  average  cloudy  skies  the  entire  tropo- 
sphere loses  heat  through  radiation.  This  is  because  of  the  nature  of  the  vertical  distribution  of  water  vapor 
in  the  atmosphere,  and  the  fact  that  water  vapor  and  liquid  water  absorb  and  emit  radiational  energy  much 
more  effectively  in  the  infrared  and  near  infrared  than  in  tbe  visible  part  of  the  spectrum.  For  all  latitudes, 
the  heat  loss  is  a maximum  in  the  middle  troposphere  for  both  clear  and  cloudy  slues. 


6.1.  CLEAR  SKIES 

The  level  of  maximum  net  cooling  in  the  case  of  clear  skies  varies  from  about  10  km  near  the  equator  to 
about  3 km  at  latitudes  60°  to  70°  N (Fig.  18).  The  value  of  this  maximum  cooling  decreases  from  about 
1.5°C  day-1  near  the  equator  to  about  0.8°C  day-1  in  the  more  northerly  regions.  In  the  lowest  2 km,  the 
rate  of  cooling  at  most  latitudes  is  quite  constant  at  about  0.7°C  day-1.  The  lower  layers  in  the  subtropical 
region  cool  at  a slightly  higher  rate  (about  0.9°C  day-1)  because  the  rapid  decrease  of  water  vapor  with  height 
permits  more  infrared  radiation  to  be  lost  to  space.  Aloft,  near  the  tropopause,  the  net  cooling  varies  from 
about  1.0°C  day-1  (0.10°  N)  to  about  0.5eC  day-1  (60°  to  70°  N).  Thus,  contrary  to  the  results  of  Moller 
(1935)  but  in  support  of  those  of  Albrecht  (1950)  it  is  found  that,  even  in  the  case  of  dear  skies,  radiational 
processes  lead  to  the  development  of  thermal  instability  in  the  lower  layers  of  tbe  tropoephere. 

6X  AVERAGE  CLOUDINESS 

The  pattern  of  net  cooling  is  somewhat  modified  if  the  distribution  of  average  cloudiness  is  considered 
(Fig.  19).  The  level  of  maximum  cooling  is  lowered  at  all  latitudes  except  near  the  equator.  Here,  above 
about  3 km,  the  air  cools  at  a fairly  constant  rate  of  about  1.2°  to  1.4°C  day-1.  The  approximate  level  of 
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Fig.  19.  Net  radiational 
cooling  in  the  tropoephere 
(°C  day'1),  average  cloudi- 
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maximum  cooling,  which  is  indicated  by  the  heavy  dashed  line  in  Fig.  19,  varies  from  about  5 km  at  10s  to 
20°  N to  about  2.5  km  at  60°  to  70°  N.  The  rate  of  net  cooling  at  this  level  is  about  1.8°C  day-1. 

The  most  striking  and  climatologically  important  difference  between  the  net  cooling  distribution  for 
clear  skies  and  that  for  a troposphere  containing  average  cloudiness  is  that  the  cooling  at  the  middle  cloud 
level  has  almost  everywhere  been  increased  — in  the  polar  regions  this  cooling  has  actually  been  doubled 
as  a result  of  the  presence  of  clouds  — and  the  cooling  in  the  lowest  kilometer  has  been  markedly  decreased. 
This  produces  a gradient  of  temperature  change  which  is  about  two  or  three  times  larger  for  average  cloudi- 
ness than  for  clear  skies.  The  marked  effect  of  cloudiness,  therefore,  is  to  intensify  the  development  of 
radiational  instability  in  the  layers  below  the  level  of  the  middle  clouds. 

Above  the  level  of  maximum  cooling,  the  gradient  of  temperature  change  is  of  opposite  sign  and  weaker 
than  the  gradient  below  this  level.  This  results  in  a stabilization  of  the  upper  levels  and  is  independent  of 
cloudiness  in  the  troposphere. 


6.3.  THE  AVERAGE  RADIATIONAL  COOLING  OF  THE  TROPOSPHERE 

The  average  radiational  temperature  change  in  the  troposphere  can  be  computed  by  considering  the 
radiation  flux  divergence  between  the  surface  and  the  tropopause  for  each  latitude  belt.  The  distribution 
with  latitude  of  this  average  radiational  temperature  change  is  shown  in  Table  12.  Columns  a and  b repre- 
sent the  cooling  due  to  infrared  radiation  and  the  heating  due  to  absorption  of  insolation,  respectively.  A 
comparison  of  values  in  Table  12  for  clear  and  cloudy  skies  shows  that  the  average  infrared  cooling  is  not 
affected  by  cloudiness,  except  at  higher  latitudes  where  the  infrared  cooling  has  been  increased  slightly  as  a 
result  of  the  clouds  being  concentrated  at  levels  below  3 km.  The  heating  due  to  absorption  of  insolation 
is  not  changed  at  all  as  a result  of  cloudiness.  This  results  from  the  fact  that  the  transmission  of  insolation 
is  decreased,  although  there  is  increased  absorption  at  the  cloud  levels;  thus  the  absorption  below  these  cloud 
levels  is  diminished.  The  net  result  is  an  average  cooling,  approximately  constant  with  latitude,  of  about 
1°C  day-1.* 

During  the  month  of  March  the  average  temperature  in  the  troposphere  is  rising,  albeit  slowly.  The 
energy  required  to  raise  the  temperature  by  the  observed  amount  is  very  small  when  compared  to  the  energy 
lost  by  the  atmosphere  through  radiation  (Panofsky  (1950)).  Compensating  processes  must  be  present, 
therefore,  to  balance  the  radiational  heat  loss.  These  processes  are  sensible  heat  transfer  through  eddy 
conduction  from  the  earth's  surface  and  the  release  of  latent  heat  by  condensation. 

In  a stable  atmosphere,  the  vertical  transport  of  heat  due  to  turbulence  is  small,  and  is  directed  down- 
ward (e.g.,  Haurwitz  (1941),  Robinson  (1951)).  Surface  heating  during  the  daytime,  however,  develops 
instability  in  the  lower  layers  which  aids  the  turbulent  exchange  of  heat.  Since  the  air  is  then  unstable, 
the  transport  of  heat  will  be  upward  to  the  top  of  the  convective  layer.  Thus  the  vertical  transport  of 
sensible  heat  is  aided  when  that  transport  is  directed  upward  and  suppressed  when  that  transport  is 
directed  downward. 

The  source  of  atmosphere  water  vapor  is  at  the  earth’s  surface  (the  oceans).  The  vertical 
transport  of  latent  heat  is  therefore  always  directed  upward  and  this  heat  is  realized  by  the  atmos- 
phere through  condensation. 


* It  is  interesting  to  note  that  this  value  of  the  heat  loss  by  radiation  in  the  troposphere  is  the  same  as  the  value  computed  by 
Gold  about  40  years  ago.  (See  the  discussion  on  Radiation  in  the  Centenary  Proceeding gs  of  the  Royal  Meteorological  Society,  1950.) 
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-1.4 
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-1.6 
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-1.44 
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-1.0 

-1.4 

+0.4 

-1.0 

40*-40* 

-l.S 

+0.44 

-0.44 
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+0.48 

-0.08 

40*-80* 

-1.1 

+0.4 

-0.4 

-1.4 

+0.8 

—0.9 

w-ao* 

-0.9 

+0.S 

-0.7 

-1.1 

+0.4 

-0.9 

40M0» 

-0.44 

+0.1 

-0.44 

-0.0 

+0.1 

-0.8 

It  Is  of  eome  intaraat  to  oompare  eetimate*  of  magnitude  of  tha  energy  involved  for  the  heat  aouroea  and 
ainlu  for  tha  atmoaphara.  Tha  net  radiational  flux  divergence  oorreaponding  to  average  cloudy  aide*  during 
March  for  a column  of  unit  croaa  aaetfoo  in  the  tropoaphera  variaa  from  288  cal  cm”*  day”'  at  0*  to  10*  N 
(depth  16  km)  to  158  cal  cm"1  day*'  at  60*  to  70*  N (depth  10  km).  If  thie  flux  divergence  ia  weighted  by 
tha  relative  area  of  each  latitude  bait,  tha  average  radiational  heat  loaa  during  March  for  the  entire  tropo- 
aphare  (0*  to  70*  N)  ia  204  cal  cm*1  day*1. 

The  quantity  of  heat  wlcaaed  to  tha  atmoaphara  by  oondonaatkm  can  be  eetimated  from  exiating  pub* 
Halted  data.  If  It  la  aaaumad  that  there  la  no  change  in  average  oloudineae  amount  through  tha  month  of 
March,  the  total  quantity  of  water  vapor  which  ia  ooodcnacd  in  tha  atmoephere  folia  to  tha  ground  aa  precipi- 
tation. Than  tha  obeerved  diatribution  of  prodpltatioa  with  latitude  ia  a meaaure  of  tha  heat  reieaaed  by 
oondenaation,  aaauming  that  there  ia  no  divergence  of  tha  horieontal  flux  of  liquid  water  acroea  each  lati- 
tude bah. 

Tha  aonal  diatribution  of  annua)  predpltation  ia  given  by  Raurwitx  and  Auatin  (1944)  after  tha  data  of 
Meinardua.  For  tha  region  0*  to  70*  N,  the  annual  amount  of  precipitation  ia  given  aa  about  105  cm,  or 
about  8.8  cm  for  March  if  that  month  ia  regarded  aa  average  with  reepect  to  precipitation.  Haurwita  and 
Auatin  point  out  that  tha  amount  of  precipitation  taken  from  Meinardua'  data  ia  probably  a little  too  high, 
ainoe  tha  precipitation  foiling  over  tha  ocean  areae  haa  bean  overea  time  ted.  Jacob*  (1951)  givea  value*  of 
tha  aonal  diatribution  of  precipitation  over  tha  oceana  for  each  aeaaon.  Hia  value*  indicate  that  a more 
oorreot  amount  of  tha  predpltation  foiling  in  tha  region  0*  to  70*  N during  March  la  about  7.5  cm.  Thia 
oorreeponda  to  an  energy  relearn  to  tha  atmoaphara  of  about  145  cal  cm*1  day”1,  if  all  of  the  condenaation 
oooura  at  a temperature  of  —10*C.  Thu*,  heat  of  oondenaatioa  offteta  about  70  percent  of  the  radiative 
boat  loaa. 

There  are  very  few  reliable  data  on  tha  average  eddy  flux  of  aenaibie  beat  from  the  ground  to  the  atmoe- 
phere.  From  tha  vahtee  calculated  for  tha  radiational  energy  loaa  and  that  eatimated  for  the  latent  heat 
gained  by  tha  atmoaphara,  it  ia  aeen  that  tha  conduction  of  heat  to  tha  atmoaphara  thould  be  of  the  order  of 
about  60  oal  cm**  day*1.  Thia  would  than  repreeent  tha  remaining  80  percent  of  the  energy  naoaaaary  to 
balance  tha  energy  loot  through  radiation. 

U.  CONCLUSIONS 

The  raauha  eta  tad  ia  tha  preceding  aaction  indicate  oiaarly  that  tha  tropoaphera  oonatitutea  an  energy 
aouroa  inaofor  aa  ha  radiative  propartiaa  are  concerned.  Although  the  net  radiative  heat  loaa  ia  fairly  uniform 


from  latitude  to  latitude,  it  is  found  that  at  all  latitudes  the  largest  radiative  cooling  is  in  the  middle  tropo- 
sphere. This  is  because  of  the  major  importance  of  water  vapor  in  atmospheric  radiative  heat  transfer  and 
the  characteristic  distributions  of  water  vapor  with  height  in  the  lower  atmosphere. 

The  presence  of  clouds  does  not  qualitatively  alter  the  cooling  distribution  but  more  sharply  defines  the 
layer  of  maximum  cooling  and  increases  the  heat  loss  in  this  layer.  In  this  way  cloudiness  contributes  to 
the  development  of  thermal  instability  below  about  3 to  5 km  and  thus  aids  in  the  transport  of  heat  upward 
through  eddy  conduction  of  sensible  heat  and,  more  important,  through  eddy  transfer  of  latent  heat. 

It  was  stated  at  the  beginning  of  this  paper  that  the  primary  function  of  the  general  circulation  is  to 
maintain  an  overall  heat  energy  balance  of  the  atmosphere.  In  the  equatorial  latitudes  the  energy  of  sensi- 
ble and  latent  heat  transported  to  the  atmosphere  is  larger  than  that  lost  by  the  atmosphere  through  radiative 
processes.  Indeed  if  one  considers  the  balance  between  total  incoming  and  total  outgoing  radiative  currents 
at  the  top  of  the  atmosphere,  there  is  a net  radiative  excess  for  all  latitude  belts  south  of  about  40°  N (Baur 
and  Philipps  (1935),  London  (1951)). 

In  the  polar  regions,  where  the  vertical  transport  of  both  sensible  and  latent  heat  is  smaller  than  the 
radiative  heat  loss,  it  is  found  that  there  is  a large  radiation  deficit  between  the  total  upward  and  total  down- 
ward radiative  currents  at  the  top  of  the  atmosphere.  Thus  the  mechanism  which  describes  the  maintenance 
of  an  equilibrium  heat  budget  for  the  troposphere  must  account  for  the  predominantly  vertical  transport  of 
heat  in  the  equatorial  and  subtropical  regions  and  horizontal  (northward)  transport  in  mid  and  polar  latitudes. 

The  atmospheric  circulation  patterns  which  would  promote  such  a redistribution  of  heat  energy  are 
different  in  equatorial  than  in  midlatitude  regions.  In  the  equatorial  and  subtropical  regions,  the  heat  trans- 
port presumably  takes  place  through  the  existence  of  large  semipermanent  meridional  cells  as  first  described 
by  Hadley  (1735).  In  midlatitudes,  however,  the  poleward  transport  of  heat  energy  is  probably  accom- 
plished through  the  large-scale  eddies  which  are  associated  with  the  major  cyclones  and  anticyclones  in  the 
atmosphere  (Jeffreys  (1926)).  Any  proposed  energy  budget  for  the  atmosphere  should  at  least  be  qualita- 
tively consistent  with  these  circulation  patterns. 
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